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CHAPTER  I 

INTRODUCTION 

1.1.  Purpose 
The  properties  of  ferromagnetic  solids  have  interested 
physicists  for  many  years,  and  are  of  sufficient  complexity  that  a 
comprehensive  theoretical  treatment  is  not  yet  a  reality.  The  pur- 
pose of  this  dissertation  is  to  apply  the  unrestricted  (or  spin- 
polarized)  Hartree-Fock  method  to  such  a  solid.  This  scheme  has  had 
a  certain  amount  of  success  in  the  explanation  of  atomic  properties, 
and,  as  v±n.  be  shown  here,  appears  to  be  an  acceptable  model  for  a 
f erromagnet . 

Experimental  information  on  the  metallic  properties  of 
f erromagnet s,  which  are  dependent  on  the  electronic  configuration, 
has  recently  become  available,  aiid  with  it  some  analyses  to  explain 
these  properties  on  the  basis  of  the  energy  band  model.  Of  the 
ferromagnetic  materials,  nickel  is  perhaps  the  best  understood. 
Energy  band,  models  have  been  presented, "^'^  that  explain  empirically 
the  available  data  on  the  electronic  and  optical  properties  of 
nickel.  Also  within  the  last  five  years,  several  papers  have  appear- 
ed in  the  literature, 5-1°  ^ich  give  energy  band  structures  calcu- 
lated from  basic  considerations,  i.e.,  they  present  solutions  for 
Schrfldinger's  equation  in  a  crystal  using  some  form  of  one-electron 


potential.  This  potential,  in  all  cases,  has  been  derived  from 
an  atomic  calculation,  corrected  for  the  effects  of  placing  the 
atoms  on  a  crystalline  lattice. 

Although  these  previous  calculations  are  qualitatively 
quite  similar,  the  iirbitrariness  of  the  potential  used  is  enough 
to  create  differences  between  them  which  are  large  with  respect 
to  the  experimental  effects  which  are  to  be  explained.  The  calcula- 
tion described  in  this  dissertation  attempts  to  eliminate  this 
dependence  on  an  arbitrary  potential  by  solving  the  equations  self- 
consistently,  in  the  same  way  as  the  Hartree-Fock  method  used  in 
atomic  calctilations.  In  this  way,  it  is  possible  to  examine  the 
validity  of  the  approximations  which  must  be  made  in  order  to  solve 
the  equations. 

The  case  of  ferromagnetic  nickel  turns  out  to  be  extremely 
sensitive  to  slight  changes  in  these  approximations.  In  particular, 
the  form  in  vdiich  the  exchange  effects  responsible  for  the  ferro- 
magnetic structure  are  inserted  into  the  theory  can  radically  change 
the  final  results.  This  effect  was  not  pointed  out  in  previous 
calculations  on  nickel  or  other  materials,  either  because  the  solu- 
tions were  not  carried  to  self-consistency  or  because  the  particu- 
lar case  was  not  sensitive  enough  to  show  a  definite  discrepancy 
with  experiment.  In  a  sense,  then,  a  calculation  on  nickel  is  a 
test  case,  vrtiich  examines  the  validity  of  an  approximation  >diich  has 
been  used  in  many  calciilations  of  this  type. 

The  procedures  used  in  this  work  are  outlined  in  the 
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following  section,  and  described  in  more  detail  in  Chapter  II.  The 
results  and  comparison  with  the  experimental  data  are  presented  in 
Chapter  III  followed  by  the  conclusions  and^ discussion  in  Chapter  IV. 

1.2.  Oi'itline  of  Calculation 

The  one-electron  model  for  a  crystal  consists  of  the  asstunp- 

tion  that  the  electronic  wave  functions  ^  satisfy  a  Schrodinger  equa- 
tion of  the  form 

where  V(r)  is  a  potential  identical  for  all  electrons,  and  includes; 

(i)  an  attractive  term  due  to  the  nuclei  situated  on  the 
lattice  sites, 

(ii)  a  repulsive  term  due  to  the  Coulomb  interaction  with  the 
other  electrons,  and 

(iii)  an  attractive  term  which  similutates  the  exchange  effects 
of  other  electrons. 

As  a  first  approximation,  V(r)  can  be  generated  from  the 

assumption  that  the  electronic  wave  functions  are  unchanged  from  their 

atomic  values.  The  atomic  potentials  on  the  appropriate  crystalline 

lattice  sites  are  then  overlapped  to  form  what  has  been  called  a 

"superimposed-atom  potential".  This,  of  course,  makes  V(r)  dependent 

on  the  pau-ticular  atomic  configuration  chosen.  Although  the  most 

logical  choice  would  be  the  groirnd  state  configuration  of  the  atom, 

this  is  not  always  the  best  choice.  It  is  true  that  the  wave  functions 

do  not  change  greatly  on  going  from  a  free  atom  to  a  crystalline 

environment,  but  the  effective  occupation  number  of  each  type  of 

orbital  (s,p,d,  etc.)  may  change.  The  reason  for  this  is  that  the 
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discrete  atomic  levels  which  the  electrons  occupy  in  a  free  atom  broad- 
en into  bands  when  these  atoms  come  together  to  form  a  solid.  If  these 
bands  happen  to  overlap,  as  they  do  in  many  metals,  then  electrons  in 
one  band  will  "spill  over"  into  the  other,  thus  reducing  the  effective 
nvmiber  of  the  first  type  in  favor  of  the  second.  This  occurs  in  the 
case  of  nickel  (figure  1.1.),  in  which  the  overlap  of  the  3d  and  4s 
bands  alters  the  configuration  such  that  the  effective  nvmiber  of 
"d-like"  electrons  is  increased  from  the  atomic  ground  state  value  of 
8,  to  approximately  9  in  the  solid. 

However,  in  a  self -consistent  calculation,  vrtiere  the  potential 
is  regenerated  after  each  iteration,  the  arbitrariness  in  the  choice 
of  V(r)  disappears.  The  superimposed-atom  potential  is  used  only  as 
a  starting  point  for  the  self-consistent  procedure. 

Once  the  initial  potential  is  chosen,  the  energy  bands  are 
calculated  for  a  selected  number  of  vectors  in  reciprocal  space  by 
means  of  the  augmented  plane  wave  (APW)  method. -^"-^  This  method  has 
been_suff iciently  developed  and  checked  out  against  other  methods  -^     ° 
so  that  it  gives  a  solution  (for  a  given  potential)  as  accurate  as 
desired.  In  its  simplest  form,  the  AF¥  method  solves  equation  (l.l) 
for  a  V(?)  of  the  "muffin-tin"  type,  i.e.,  spherically  symmetric  in 
spheres  about  each  lattice  sitenand  constant  between  these  spheres. 
This  is  not  a  necessary  restriction,   '  but  it  can  be  removed  only  at 
cost  of  much  greater  complexity,  without  much  improvement  in  the  final 
results. 

The  energy  eigenvalues  obtained  in  this  way  are  then  used  to 
find  the  associated  wave  functions.  The  total  charge  density  is 
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generated  from  these  wave  functions,  by  summing  over  the  occupied 

levels.  The  procedxire  for  forming  the  new  muffin-tin  potential  to 

17  18 
be  used  for  the  next  iteration  is  straightforward,      with  a 

certain  amount  of  care  necessary  in  the  choice  of  the  intersphere 

constant.  The  iterations  are  continued,  going  through  the  above 

manipulations  each  time,  until  there  is  no  significant  difference 

between  the  energies  obtained  from  two  sucessive  iterations. 

Each  of  these  steps  is  described  in  detail  in  Chapter  II, 

and  the  entire  Ccilculation  represented  by  the  flow  chart  in  figure 

1.2. 
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Figure  1.2.  Flow-chart  for  the  self -consistent  calculation 


CHAPTER  II 

METHOD  OF  CALCULATION 

2.1.  The  Unrestricted  Hartree-Fock  Method 
The  Hartree-Fock  formalism  consists  of  approximating  the 
wave  fimction  for  an  N-electron  system  by  a  single  antisyiiimetrized 
product  of  one  electron  fxmctions,   u^-  of  space  and  spin  coordinates 
Xi  =  (i^jSi) 

^^    c;>{^ZL,(x,)u^(xJ  ...   VL^iX^)}   .  (2.1) 

For  a  Hamiltonian  consisting  only  of  electrostatic  interac- 
tion and  kinetic  energy  terms,  the  variational  principle  leads  to 

19 
N  coupled  integro-differential  equations  of  the  ibrm 

[II  S/>ins]  £  t  =  ■/,  2,  .  .  •  A/J   ,      (2.2) 

i^ere  <5K  is  the  one-electron  kinetic  and  nuclear  potential  operator. 
The  second  summation  is  restricted  to  one-electron  states  (j)  with 
spin  parallel  to  state  (i),  whereas  the  first  is  over  all  the  elec- 
trons including  the  ith  state  itself. 


These  equations  can  also  be  written  in  the  form 


(2.3) 


where  V  (if)  = 


-  ^ 


is  the  nuclear  potential  term. 


(atan) 
(solid) 


is  the  Coulomb  potential,  and 

J'  1  It' 

[  spin    sj 

is  the  "exchange  potential".  Note,  that  V  and  V  are  "local"  poten- 
tials, i.e.,  they  appear  in  the  equations  as  a  multiplicative  ope- 

(s) 
rator  acting  on  u^,  but  V^^^  '  is  not.  It  is  this  peculiarity  of  the 

equations  irtiich  causes  the  greater  part  of  the  difficulty  in  their 

solution.  The  exchange  potential  can  be  expressed  as  a  local  operator, 

but  with  the  result  that  it  will  have  a  different  form  for  each 

orbital: 


[II  Sf>ir>s]  ' 
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The  equations  are  solved  by  an  indefinite  number  of  iterative 
steps.  Some  initial  set  of  Uj^'s  is  chosen,  V^,  and  V^  are  generated 
from  them,  and  the  resulting  differential  equations  are  solved  for 
a  new  set  of  Uj^'s.  The  procediire  is  then  repeated  for  this  new  set 
until  the  final  u^'s  are  the  same  as  the  initial  set  (i.e.,  "self- 
consistency"  is  achieved).  Fortunately,  for  an  initial  set  of  u^'s 
close  enough  to  the  self-consistent  solution,  this  process  is  conver- 
gent. 

Because  of  the  difficulties  involved  in  solving  the  three 
dimensional  integro-differential  equations  (2.3),  their  solution 
involves  a  major  computational  effort,  and  calculations  involving 
heavy  atoms  are  few.  In  practice,  the  Hartree-Fock  formalism  is 
often  modified  by  assuming  the  potentials  in  equations  (2.3)  are 
spherically  symmetric.  (Here  we  are  referring  to  the  case  of  a  free 
atom.  The  equations  for  a  solid  will  be  discussed  in  subsequent 
sections  of  this  chapter.)  This  allows  us  to  separate  the  vra,ve 
function  in  the  form 

where  y    is  a  spherical  harmonic  with  angular  momentum  quantum 
numbers  jt     and  fn   and   ^  (s)  is  a  spin  function  with  spin  quantum 
number,  a  =-  ^.   This  one  assimiption  (known  as  the  central  field 
approximation)  effects  enormous  simplification  in  the  equations,  in 
that  ^n  integrals  and  differential  equations  become  one-dimensional. 
In  particular  (2.3)  reduces  to 
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The  exchange  potential  V^^-Cr)  is  slightly  more  con?)licated  than  Vj,(r) 
involving  an  extra  summation  over  integrals  vrtiich  contain  Pj  (r)  as 
well  as  P^Cr),  (See  reference  19,  page  17.) 

In  general,  Pj  will  depend  on  1,  m,  and  s,  as  well  as  the 
principal  quantimi  number  n  which  labels  the  eigenvalues  of  (2.6)  in 
increasing  order.  In  the  so-called  "conventional"  Hartree-Fock 
scheme,  P^  is  assumed  independent  of  m  and  s,  restilting  in  a  degene- 
racy of  2  (2>£+l)  for  each  energy  eigenvalue,  thus  reducing  the 
number  of  equations  to  be  solved.  This  has  the  advantage  of  giving 
a  final  wave  function  (2.1)  which  is  an  eigenfunction  of  the  total 
angular  momentum  L*  and  total  spin  operator  S  .  However,  this  wave 
function  may  give  entirely  false  results  in  interpreting  results  in 
\diich  spin  effects  are  important.  For  example,  consider  a  transition 
element,  for  which  the  groimd  state  has  two  4s  electrons  and  an 
unfilled  3d  shell.  The  conventional  Hartree-Fock  solution  in  this 
case  WDiold  have  a  zero  spin  density (i.e.,  the  difference  between  the 
charge  densities  of  opposite  spins)  near  the  nucleus.  But  the  inter- 
pretation of  the  hyperfine  structiire  of  an  atom  with  non-zero  nuclear 
spin  demands  that  the  spin  density  be  non-zero.  This  sort  of  effect 
is  the  justification  for  the  use  of  "unrestricted"  Hartree-Fock 
sbheme.^*^*21  -phe  terra  "unrestricted"  refers  to  the  relaxation  of 
the  restriction  that  Pi  be  independent  of  s.  The  number  of  equations 
to  be  solved  is  doubled,  giving  a  (2-£+l)  -  fold  degeneracy  for 


.12 


each  ^x\£b,   ^°^  *  configuration  iriiich  has  an  unequal  number  of  up 
and  down  spins,  such  as  an  unfilled  3d  shell,  the  relaxation  of 
the  constraint  will  result  in  a  non-zero  spin  density.  This  desira- 
ble feature  of  the  wave  function  is  gained  at  the  sacrifice  of  the 

2 

symmetiy  requirement  that  (2.1)  be  an  eigenf unction  of  S  .  However, 

22 
it  has  been  shown   that  the  unrestricted  wave  functions  are  in 

substantially  better  agreement  with  experimental  data.  It  also  should 

be  noted  here  that  the  unrestricted  solution  can  be  made  an  eigen- 

P  23 

function  of  S  by  the  application  of  projection  operators,    but 

such  a  complication  is  beyond  the  scope  of  the  present  calculation. 
Both  the  conventional  and  the  restricted  forms  of  the  Hartree-Fock 
must  be  regarded  as  approximations  whose  error  is  impossible  to  esti- 
mate without  doing  the  actual  exact  calculations.  Since  these  are 
intractable  at  the  present  time  for  atomic  systems  with  more  than  a 
few  electrons,  appeal  must  be  made  to  experiment  to  judge  the  use- 
fulness of  the  different  forms.  In  particular,  as  will  be  seen,  the 
conventional  scheme  is  much  less  suitable  than  the  unrestricted  for 
the  case  of  a  ferromagnetic  solid, 

2.2.  Exchange  Approximations 

The  exchange  potential  of  equation  (2.4)  for  a  particiilar 

orbital  u.  involves  a  sum  over  all  electrons  u.   of  spin  parallel  to 
1  J 

that  of  vi^,  of  integrals  involving  both  u^  and  Uj.  Ixjoking  at  the 
situation  from  a  one-electron  point  of  view,  this  means  that  each 
electron  moves  in  a  different  effective  potential.  For  a  large 
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nxmber  of  electrons,  this  results  in  formidable  computational 
difficiilties  even  >Aien  the  central  field  approximation  is  used. 
These  difficulties  are  not  insurmoruitable  for  a  free  atom,  but  are 
virtually  impossible  to  overcome  for  the  case  of  a  solid.  Also,  we 
are  interested  in  solving  the  problem  within  the  framework  of  the 
one-electron  model,  in  which  all  electrons  move  in  the  same  effective 
potential.  This  simplification  has  enjoyed  a  great  deal  of  success 
despite  its  obvious  limitations,  and  represents  the  present  limit 
of  complexity  to  which  calculations  in  solids  can  be  carried. 

What  is  clearly  indicated,  then,  is  the  formation  of  some 
kind  of  average  exchange  potential  which  will  retain  the  main  fea- 
tures of  equation  (2.4)  One  such  average  has  been  suggested  by 

24 
Slater,   i.e.,  ^ ,   x, 

1/  (p)  =     z_j  y  —   . 

" '     , V  ,    y  ^f^-,)%i^,)      ,  (2.7) 

[n  Sf>ins]  J 

which  is  simply  the  sum  of  V^  ^^  weighted  by  the  probability  factor, 

J 

There  is  one  particular  case  for  which  the  expression  (2.7) 
can  be  evaluated  exactly,  viz.,  the  free-electron  gas.  If  the  free- 
electron  wave  functions  11^=  exp(iki.r)  are  substituted  into  (2,7), 
\^  the  integral  can  be  easily  performed,    and  the  sum  over  k.  up  to 

the  Fermi  level  gives 
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trfiere^  is  the  constant  electron  charge  density  (in  imits  of  the 
electronic  charge). 

In  an  atom,  of  course,  the  electrons  are  far  from  free,  and 
the  charge  density  is  not  a  constant,  but  if  we  take  the  extreme 
step  of  assuming  the  exchange  interactions  between  electrons  in  an 
atom  (  or  a  solid)  depend  only  on  the  local  electronic  charge 
density,  then  a  function  is  obtained  vrtiich  tioms  out  to  be  .remarkably 
similar  to  (2.7).  Explicitly  this  function  is 

yt^^-^-  -^/4  t'\^^.^i^f.        (2.9) 

In  at  least  one  case,  that  of  the  ion  Cu"*",  a  detailled 
comparison  of  the  approximations  (2.7)  and  (2.9)  with  the  exact 
Hartree-Fock  expression  (2.4),  has  been  worked  out.  ^^  The  results 
show  that  the  Slater  average  (2.7)  and  the  free-electron  average 
(2.9)  are  practically  idential.  At  least,  the  deviations  of  (2,7) 
from  (2.9)  are  much  less  than  the  deviations  of  either  from  the 
Hartree-Fock  value  (2.4).  The  point  to  be  drawn  here  is  that  there 
is  no  advantage  in  using  (2.7)  instead  of  (2.9),  and  that  if  one  is 
forced  by  practical  considerations  into  using  a  local  exchange  poten- 
tial, then  the  errors  involved  can  be  ascribed  to  the  averagijig  pro- 
cedure rather  than  the  use  of  the  free-electron  approximation. 
Recently,  Kohn  and  Sham  ^^  by  the  use  of  a  different  averaging  scheme 
have  concluded  that  (2.9)  should  be  multiplied  by  a  factor  of  2/3. 
Calculations  on  various  atoms  ^7,28  ^^^^  ^^^.^  ^^^^   perhaps  this  is 
closer  to  the  truth,  in  that  the  resultant  wave  functions  are  more 
similar  to  the  Hartree-Fock  values.  Undgren  ^9,30  ^^  ^^^^  ^ 
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pcirameterized  one-electron  exchange  potential  and  minimizing  the 
correct  total  Hartree-Fock  energy  vdth  respect  to  these  parameters 
gets  amazingly  good  agreement  with  the  best  Hartree-Fock  values. 
This  "optimum"  potential,  he  finds, ^  is  always  smaller  in  magnitude 
than  (2.9).  The  consensus  appears  to  be  that,  although  the  free- 
electron  approximation  V^  reproduces  the  general  features  of  the 
Hartree-Fock  scheme,  it  overestimates  the  exchange  effects. 

The  extension  of  this  approxijnation  to  the  unrestricted 
case  is  straightforward.  In  the  coxirse  of  the  derivation  of  (2.8), 
half  of  the  electron  gas  was  assimied  to  have  up  spin,  and  the  other 
half  down  spin.  This  introduces  an  extra  factor  of  2  in  equation  (2.9) 
iriien  we  replace  the  total  charge  density  by  that  of  only  one  spin. 
The  effective  exchange  potential  is  then 

where  the  svunmation  now  contains  the  orbitals  of  one  spin  only.  This 
is  the  approximation  used  in  this  calculation.  Substitution  of  (2.10) 
into  (2.6)  gives  the  actual  equations  to  be  solved: 

where  V^U)  ^  -Zi;    ^  3:.   f  (rirO  df' -h  1 J  (rC>-0^dy^ 

^-t  (2.11) 
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Wjj*_  is  the  number  of  electrons  in  ths  conf igiiration  with  quantum 
ntimbers  n  andJl,   and  spin  s. 

These  equations  are  exactly  the  same,  vfith  the  exception 
of  the  added  spin  subscript,  as  those  used  by  Herman  and  Skillman 
in  their  atomic  calculations.    Their  program,  vrtiich  was  published 
with  their  calculations,  was  modified  for  the  unrestricted  case 
during  the  course  of  this  work.  As  a  test  case,  this  program  was 
used  to  find  the  wave  function  for  the  Fe  atom  in  the  configuration 

(3do()^  (3dj5)"^  (ifSd)"^  (4sB)  ,  in  order  to  compare  with  the  results 

21 
of  Wood  and  Pratt,    who  did  exactly  this  calculation.  A  comparison 

of  the  one-electron  energies  is  shown  in  table  2.1.  The  energy  values 
in  the  unrestricted  case  tend  to  bracket  those  of  Herman  and  Skillman 
and  the  splittings  are  larger  for  the  rF=3  electrons  than  for  the 
other  shells,  as  wotad  be  expected.  Apart  from  a  large  discrepancy 
(almost  certainly  an  error  in  their  paper)  at  the  Is  level,  the 
results  are  comparable  to  those  of  Wood  and  Pratt.  The  results  are 
not  identical  in  that  their  calculation  was  not  taken  out  to  as 
great  a  degree  of  self-consistency. 

At  first  glance,  the  discrepancy  between  the  energies  calcu- 
lated using  the  free  electron  approximation  and  the  exact  Hartree- 
Fock  values  seems  discouragingly  large.  However,  this  does  not  mean 
that  the  corresponding  wave  function  is  inaccurate  to  the  same  degree. 
As  Lindgren^  has  pointed  out,  Koopman's  theorem  no  longer  holds 
when  the  Hartree-Fock  equations  are  solved  approximately,  and 
therefore  the  one-electron  energies  no  longer  have  the  same 
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Table  2.1.  Hartree-Fock  one  electron  energy  parameters  for 
Fe  (3cic()5  (3d^)3  (4s)2 


Exact 
CHF^ 

Approx. 
CHF^ 

Approx. 
UHF<^ 

Approx. 
UHP^ 

^ 

fi 

o( 

3 

Is 

-522.7® 

-515.8 

-584.5 

-584.2 

-515.8 

-515.8 

2s 

-63.84 

-60.96 

-61.08 

-60.40 

-61.03 

-60.75 

2p 

-54.79 

-53.08 

-53.16 

-52.65 

-53.12 

-52.91 

3s 

-8.308 

-7.269 

-7.463 

-6.930 

-7.492 

-6.929 

3p 

-5.455 

-4.891 

-5.061 

-4.540 

-5.111 

-4.558 

3d 

-1.271 

-0.963 

-1.122 

-0.664 

-1.168 

-0.664 

4s 

-0.510 

-0.545 

-0.532 

-0.428 

-0.592 

-0.491 

^  Exact  conventional  Hartree-Fock  values,  taken  from  Watson^^, 

b  "^^ 

Approximate  conventional  Hartree-Fock  values  (V^^  —  ^x  ^  >   taken 
from  Herman  and  Skillman32, 

^   Approximate  unrestricted  Hartree-Fock  values  (V^^  =  V*^)  ,  taken 
from  Wood  and  Pratt21, 

^  Approximate  unrestricted  Hartree-Fock  values  (V^^  =  V**)  ,  calcu- 
lated during  the  course  of  this  work. 

®  All  energies  are  in  Rydbergs. 


18 


significance  of  being  equal  to  the  binding  energies,  ^y-  applying 
the  appropriate  corrections,  he  has  shown  that  binding  energies 
can  be  extracted  from  the  approximate  solutions  which  are  in  much 
better  agreement  with  the  exact  values. 

To  further  examine  the  method,  a  calculation  was  done  for 
the  ion  Ni"^  in  the  configuration  Odd)^   (3d/)3  ,  for  vrtiich  an  exact 
unrestricted  Hartree-Fock  solution  is  available.    The  results  are 
shown  in  table  2. 2, both  for  the  standard  exchange  and  for  one  2/3 
as  large.  Note  that  the  discrepancies  with  the  exact  energy  values 
are  larger  in  the  latter  case,  but  that  core  polarisation  (figure  2.1) 
is  in  better  agreement.  This  serves  as  an  illustration  of  the  fact 
that  inaccuracies  in  the  one-electron  energies  are  not  necessarily- 
indicative  of  inaccuracies  in  the  wave  function  itself. 

2.3.  Superlmposed-Atom  Potential 
In  the  preceding  section  we  have  shown  how  to  obtain  a 
one-electron  potential  valid  for  electrons  in  an  atom.  We  are  now 
faced  with  the  problem  of  how  to  extend  this  method  when  the  atoms 
come  together  to  form  a  solid.  Suppose  we  make  the  rash  assumption 
that  the  electron  charge  density  in  the  solid  is  unchanged  from  the 
atomic  values.  If  />,  (r)  is  the  spherical  charge  density  of  spin  s 
in  the  atom,  derived  from  a  self -consistent  calculation,  then  the 
charge  density  in  the  solid  i^(r)  is  just  a  sum. 


/,(?;  =  j>j^)  .  ^y'^  (/"-?,'}. 


(2.15) 
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Table  2.2.  Hartree-Fock  one-electron  energy  parameters  for 
Ni+^  (3do()5  (3d^)3 


Exact 

UHF^ 

Vx 

Approximate 
^x 

UHF^ 
Vx  = 

3 

-612.6  c 

3                         o( 

j3 

o( 

^ 

Is 

-612.6    -605.3 

-605.3 

-597.2 

-597:2 

2s 

-77.30 

-77.15    -74.20 

-74.04 

-72.17 

-72.08 

2p 

-67.35 

-67.21    -65.47 

-65.35 

-63.27 

-63.20 

3s 

-11.35 

-10.96    -10.19 

-9.901 

-9.477 

-9.294 

3p 

-8.200 

-7.725    -7.445 

-7.158 

-6.762 

-6.583 

3d 

-2.898 

-2.705    -2.782 

-2.519 

-2.266 

-2.m 

a 

Exact  unrestricted  Hartree-Fock  values,  taken  from  Watson  and 
Freeman,20 

^  Approximate  unrestricted  Hartree-Fock  values,  calculated  during 
the  course  of  this  work, 

°  All  energies  are  in  Rydbergs. 
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where  a^^  is  the  position  vector  of  the  ith  lattice  site.  We  consider 

a  particular  case  here,  that  of  only  one  atom  in  a  imit  cell 

The  Coulomb  potentials  defined  in  (2.11)  will  also  combine  linearallyj 


1/ 


(j?}  =  \/^{H    +  2^    Vcd^'^iO 


(2.13) 


where  Vc(r)  is  the  atomic  Coulomb  potential  and  the  exchange  is 

\si^^-    -   <^[i-^^s(^>f    ■  (2.14) 

For  computational  reasons,  vrtiich  will  be  made  clear  in  the 
next  section,  we  wish  to  have  a  potential  which  will  have  the  so- 
called  muffin-tin  form,  i.e., spherically  symmetric  in  non-overlapping 
spheres  of  radius  R  about  each  lattice  site,  and  a  constant  between 
these  spheres.  To  find  the  spherical  average  ^J^^    of  the  charge 
density  (2.12),  we  integrate  over  all  angles  to  get 

I  0 

<»-i  +  >'  (2.15) 


-j'^^'-^'-Z^J''^'''^^'' 


I  Art- 1 


where  nj^  is  the  number  of  neighboring  atoms  at  distance  a^. 

The  spherically  averaged  Coulomb  potenticil^V^  has  the 
same  form,  and  by  substituting  (2.15)  into  (2.14)  we  can  approximate 
the  spherically  averaged  exchange  potential  ^V^^^  .  _  . 

The  constant  intersphere  potential  V   is  defined  by 


(2.16) 
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vrtiere  SI    =  J^  —  kTT    R   is  that  part  of  the  unit  cell  SI  which 
lies  outside  the  spheres.  This  can  be  ev£j.uated  in  terms  of  the 
spherically  averaged  quantities  by  splitting  up  the  integral  over 
Vq  (r)  into  separate  parts  overji  and  the  sphere 72^,  i.e., 

.00  f.^  _^ 


=    4-Tr  ff^Vjt)  oir    -    47r|A'»<l/> 

a  U 


0 

finally  we  have 


/t 


^«^«      Qo  =  Ji^  '  "^^  [ ''' <fs>  ^'^ 

is  the  charge  of  spin  s  in  the  volume  J2^  and  ^   is  the  fraction 
of  electrons  with  spin  s. 

This  scheme  affords  a  straightforward  method  of  deriving  a 
one-electron  potential  suitable  for  a  solid  from  the  results  of  a 
self-consistent  atomic  calculation.  The  sum  over  neighbors  in  equa- 
tion (2.15)  typically  converges  after  two  or  three  terms,  and  the 
integrals  involved,  being  all  one-dimensional,  present  no  difficul- 
ties. The  sphere  radius  R  is  arbitrary,  but  better  convergence  in  the 
energy  values  is  obtained  in  the  solution  of  the  energy  band  problem 
if  R  is  taken  as  large  as  possible,  i.e.,  equal  to  one-half  the 
nearest  neighbor  distance.  For  the  face  centered  cubic  structure  with 
one  atom  per  unit  cell,  ft=  a/2V?,  where  a  is  the  "lattice  parameter". 
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the  distance  between  second  nearest  neighbors. 

There  are  two  basic  difficulties  with  the  potential  descri- 
bed here: 

(1)  The  neglect  of  the  non-sphericeil  terms  in  the  potential. 
This  can  be  remedied  by  adding  them  in  as  a  perturbation  after  a 
solution  is  found.  However,  if  these  terms  are  not  negligible,  formi- 
dable complexities  are  introduced  into  the  solution. 

(2)  The  assumption  that  the  charge  density  is  unchanged  from 
its  atomic  value  is  not  always  valid.  The  orbitals  will  certainly  be 
distorted  to  some  degree  by  the  influence  of  neighboring  atoms.  This 
would  be  much  more  noticeable  in  a  metal,  vrtiere  the  wave  functions 
have  a  substantial  overlap,  than  in  an  insulator  where  the  charge 
density  is  more  confined.  It  mav  be  possible  to  choose  an  atomic 
configuration  other  than  the  ground  state  vriiich  will  'generate  a 
charge  density  close  to  the  actual  charge  density  in  the  solid,  but 
this  is  impossible  to  determine  non-empirically. 

However,  this  problem  will  not  conceni  us  in  this  work, 
since  we  will  be  finding  a  self -consistent  solution.  The  superimposed- 
atom  potential  as  described  here  is  used  only  as  a  logical  starting  point 
for  the  first  iteration. 

2.4.  The  Energy  Band  Calculation 
Within  the  framework  of  the  one-electron  model  for  a  solid 
containing  N  atoms  and  ZN  electrons,  the  total  electronic  wave  func- 
tion is  an  ZN  by  ZN  determinant  of  functions  t*!,  each  of  which 
satisfy  an  equation  of  the  form 
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[-V^^  V(F^)]ll5,  .  E,X|5,  , 


(2.19) 


vAiere  E^  are  the  allowed  energy  levels.  The  potential  V(?)  includes 
the  interactions  vdth  the  N  nuclei,  vrtiich  are  assumed  to  be  in  fixed 
positions  on  the  lattice  sites,  and  simulates  the  interactions  vrLth 
all  the  other  electrons.  Effectively,  the  problem  of  solving  the 
Schrc)dinger  equation  for  the  (Z+l)N  particles  has  been  replaced  by 
an  easier  one,  in  vdiich  the  ZN  electrons  are  pictured  as  non-inter- 
acting particles  moving  in  a  fixed  potential.  In  this  picture,  the 
Fermi  statistics  hold,  so  that  the  ground  state  of  the  system  is 
obtained  by  choosing  the  solutions  of  (2.19)  corresponding  to  the 
lowest  ZN  values  of  E^. 

For  a  solid  with  a  ferromagnetic  structiire,  the  problem  is 
generalized  slightly,  in  that  the  electrons  are  divided  ±n  two 
separate  systems,  according  to  their  spin,  ;^*iich  interact  differently 
with  their  environment.  In  the  one-electron  model,  this  means  solving 
two  equations: 

•■  •-      '■'•''       'J  •      (2.20) 

The  lowest  ZN  eigenvalues  of  the  combined  set  of  solutions 
will  represent  the  groimd  state.  Since  the  potentifids  V.  and  Vg 
are  unequal,  the  two  sets  E^  and  E^  will  be  different.  This  will,  in 
general,  imply  an  excess  of  one  spin  over  the  other  in  the  ground 
state,  vAiich  is  found  experimentalljr  in  a  ferromagnetic  material. 
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Once  this  model  has  been  chosen,  vre  are  left  vrith  the  pro- 
blem of  solving  the  equations  (2.20).,  This  can  be  made  feasible  by 
the  application  of  the  symmetry  properties  of  the  Hamiltonian  3t  = 
-vVvgCr). 

First,  the  translational  symmetry  of  ^  allows  us  to  choose 
the  eiganfunctions  *»lJ-(r)  to  be  Bloch  functions,  i.e.,  to  satisfy  the 
relation 

-^■(^■^a.)  -  e"  '"^^   llJ^CT^)  ,         (2.21) 

where  a^  is  a  lattice  vector,  and  k  is  a  real  vector  which  takes  on 
N  values  corresponding  to  the  N  one-dimensional  representations  of 
the  translation  group.  In  the  language  of  group  theory,  a  Bloch 
function  is  a  basis  function  for  the  k-th  representation  of  this 
translation  group  T. 

The  reciprocal  lattice  vectors  g^  are  defined  by  the  rela- 
tion 


|..  Sj  =  Zffny  , 


(2.22) 


where  the  n^A   are  integers. 

If  I  k>  is  a  Bloch  function  belonging  the  k-th  representation 
of  T,  then  it  follows  from  (2.21)  and  (2,22)  that  |  k  +  g^y  also 
belongs  to  the  same  representation,  and  further  ,  that  1  k^  and  Ik''^ 
belong  to  the  same  representation  if  and  only  if  k-k'  is  a  reciprocal 
lattice  vector.  We  define  the  Brillouin  zone  as  that  set  of  N  vectors 
ko  in  k-space  satisfying  the  relation: 
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T*    -> 


k,-t.     ^  TT    y  dj  .  (2.23) 

Every  vector  in  k-space  is  equal  to  the  sum  of  a  "reduced"  vector 
ICq  (in  the  Brillouin  zone)  and  a  reciprocal  lattice  vector. 

Suppose,  we  have  available  a  complete  set  of  Bloch  functions 
I  k^  corresponding  to  all  possible  vectors  in  k-space,  then  "^^   can 
be  expressed  as  a  linear  combination  of  these  functions.  But,  if  \[r. 
is  to  satisfy  (2.21),  then  only  those  functions  which  belong  to  the 

*^  4.1. 

k   representation  will  contribute  to  it.  Therefore  iji.  can  be  expres- 
sed as  a  linear  combination  of  that  subset  of  the  complete  set  which 
has  the  same  reduced  vector,  i.e., 

is  the  most  general  eigenfunction  corresponding  to  the  irreducible 
representation  associated  with  the  reduced  vector  Ic^ . 

A  secular  equation  for  a  particular  1^  can  now  be  set  up  by 
substituting  (2.24)  into  (2.20): 

T  <?.+  f  I  3€-  £{X.)\K^f>%f    =0    (2.25) 

and  the  energies  E(kQ)  can  be  found  by  the  solution  of  the  determi- 
nantal  equation 
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and     Aj^'  =  <■?,+  Jl  it,1-f>. 

The  order  of  the  determinant  in  (2.26)  will  be  equal  to  the 
nvunber  of  lattice  vectors  needed  to  give  convergent  values  for  eCICq). 
This,  of  course,  will  be  dependent  on  the  basis  set  chosen.  For 
example,  the  plane  waves,  e-^  form  a  complete  set  of  Bloch  functions, 
but  are  quite  unsuitable  for  a  problem  in  which  the  potential  is 
atomic-like . 

For  a  "muffin-tin"  potential,  as  described  in  Section  2.3, 

suitable  basis  functions  are  the  augmented  plane  waves    (APW'S) 

defined  by 

•  "^  -^ 

I  '^^/      ~       ^  outside  the  spheres, 

inside  the  spheres, 

where  Y,  are  spherical  hannonics,  OLj^  satisfy  a  radial  equation 
corresponding  to  energy  E,  and  the  a^  are  constants  defined  so  that 
["^^is  continuous  at  the  sphere  radius  r  =  R,  (see  Appendix  I).  These 
satisfy  the  requisite  Bloch  condition  (2.21),  and  the  matrix  elements 
required  for  the  secular  equation  (2.26)  are  easily  derived.  For 
reference,  these  elements  are  listed  in  Appendix  I,  equations  (AI-l6,17). 

The  translational  symmetry  used  to  derive  the  secular  equa- 
tion (2.26)  is  only  part  of  the  total  symmetry  of  a  crystal.  In  gene- 
ral, the  Hamiltonian  is  invariant  under  a  space  group,  vrtiose  elements 
are  combinations  of  translations  and  rotations  (either  proper  or 
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improper).  The  translation  group  is  always  a  subgroup  and  in 
certain  cases,  called  synnnorphic ,  the  space  group  can  be  expressed 
as  a  semi-direct  product  of  the  translational  subgroup  and  another 
subgroup  called  the  point  group,  consisting  only  of  rotations.  This 
is  the  case  for  the  face-centered  cubic  structure,  for  vrtiich  the 
point  group  is  the  fiill  cubic  group  Oj^. 

The  representations  for  a  syimnorphic  ^pace  group  are  parti- 
cularly simple,  and  it  can  be  shown  ^^  that  a  basis  fvmction  |lf  > 
of  the  translation  group  can  also  be  chosen  to  be  a  basis  function 
for  that  subgroup  of  the  point  group  which  consists  of  elements 
transforming  the  vector  ICq  into  an  equivalent  vector.  This  subgroup, 
known  as  the  group  of  the  reduced  vector  %,   can  be  defined  mathemati- 
cally as 


^(0=   i^l^K-K^s}  , 


(2.28) 


where  ^  are  the  rotational  elements  of  the  space  group. 

A  basis  function  for  this  group  can  be  manufactured  from 
I W  ^y  ^^^   application  of  the  standard  "projection"  operator 

vrfiere  "plXjjjis  the  (l,j)  matrix  element  of  the  pth  representation 
of  J(i^.  The  projected  function. 


(?.  prj>  B    0lli> 


(2.30) 
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vdll  transform  as  the  I''"  row  of  Tp,  i.e., 

I' 

It  can  be  shovm  ^-'~-'°  that  the  effect  of  operating  on  a 
Bloch  function  I k  >  such  as  an  AFW,  with  the  rotation  3{   is  to 
transform  it  into  another  Bloch  function  corresponding  to  9^  k^ 
Therefore  the  projected  basis  function  is 

The  projection  operators  (2.29)  satisfy  the  orthogonality 
relation 

^JL  ^MN  "   ^p  ^P^   ^^  ^"^  '  (2.33) 

vrtiere  G  is  the  order  of  ^(^ )  and  n  is  the  dimension  of  "  . 
o  <l  o     p  p 

Because  of  the  relation  (2.33),  any  operator  A  which  is 
invariant  under  all  operations  ofO[(kQ)has  zero  matrix  elements  between 
projected  functions  belonging  to  different  representations,  or  to 
different  rows  of  the  same  representation.  It  follows  directly  from 
(2.33)  that 

f  <^.l  Al  K,PJ-J-'>- 

The  requisite  overlap  and  energy  integrals  are  found  by 
substituting  A=I  (the  identity  operator)  and  A=T^respectively.  The 
explicit  form  of  these  matrix  elements  between  projected  APW's  is 
shown  in  Appendix  I. 
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If  the  projected  functions  1^^  PIJ>  are  chosen  instead  of 
1 1^  ,   then  the  secular  equation  (2.26)  is  considerably  simplified 
at  symmetry  points,  i.e.,  where  (jCl^J  contains  more  than  one  element. 
The  orthogonality  of  equation  (2.34)  has  the  effect  of 
•^  "splitting"  up  the  secular  equation  into  blocks,  one  for  each  combi- 

nation  of  P  and  I,  which  can  be  solved  separately.  Also,  since  the 
right  hand  side  of  (2.34)  is  independent  of  the  row  I,  then  we  need 
solve  for  only  one  row  of  each  representation.  Thus  for  each  P  we 
have  a  smaller  secular  equation: 

Where      ^^I,.  J^^   =   <t,.^  \  ^  \  Z.^f    ^ ''^> 

The  projected  functions  (2.30)  at  a  symmetry  point  are  not 
all  independent.  In  Appendix  II  is  given  a  recipe  for  selecting  out 
of  all  the  functions  |ko  +  f  P  I  J  >,  a  linearally  independent  set. 

In  practice,  the  energies  E(ko)  are  calculated  at  a  few 
symmetry  points,  because  of  the  saving  in  tme  afforded  by  the  use 
of  proiected  functions,  and  some  kind  of  interpolation  scheme  is 

T 

"*^  used  to  find  the  energy  bands  over  the  whose  Brillouin  zone.  It  xs 

a  straightforward  matter  for  a  computer  to  calculate  the  matrix 
elements  (equations  AI-16  to  22)  and  diagonalize  the  secular  equa- 
tion. There  are  two  infinite  sums  involved  (over  the  i  values,  and 
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the  lattice  vectors  f),  but  empirically  it  is  found  "^^  that,  for 
a  metal,  6i  values  and  50  lattice  vectors  (at  a  non-symmetry  point) 
are  enough  to  give  convergence  in  the  energy  values  to  less  than 
.001  Rydberg.  This  means  that  one  never  has  to  diagonalize  a  matrix 
larger  than  50  by  50,  and  at  a  symmetry  point,  or  if  less  accxiracy 
is  desired,  considerably  less. 

2.5.  The  Self-Consistent  Calculation 
In  the  preceding  section,  it  has  been  shown  how  the  energy 
eigenvalues  and  eigenfunctions  can  be  found  from  a  given  potential, 
(  or  two  potentials  as  in  the  ferromagnetic  case).  If  a  self -consis- 
tent solution  is  desired,  there  remains  the  problem  of  generating, 
from  the  calculated  wave  functions,  a  new  potential  to  be  used  on 
the  subsequent  iteration, 

Vflien  an  approximate  exchange  of  the  type  given  in  equation 
(2.8)  isjised,  the  one-electron  potential  is  a  functional  of  the 
charge  density  p(?) 

S 

where  V^  is  the  particular  exchange  approximation  to  be  used,  and 
Vj,  is  the  Coulomb  potential  satisfying  Poisson's  equation. 


V  V^  =  -  StTJ)       (atomic  units). 


(2.37) 


In  the  computations  reported  here,  the  charge  density  is  assumed  to 
be  spherically  symmetric  in  a  sphere  about  each  lattice  site,  and 
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constant  between.  The  departwes  of  the  actual  charge  density  from 
this  simplified  model  can  be  handled  as  pertiirbations  on  the  final 
solution,  but  only  at  a  cost  of  considerable  complication. 

In  Appendix  III  it  is  shown  how  to  obtain  an  averaged  charge 
density  of  this  type  from  the  occupied  wave  functions.  The  quantities 
to  be  calculated  are  closely  related  to  those  already  found  in  the 
expression  for  the  overlap  integrals  used  in  the  secular  equation. 
Taking  advantage  of  this  fact  simplifies  the  work  to  be  done. 

The  formation  of  the  Coulomb  part  of  the  muffin-tin  poten- 
tial from  the  averaged  charge  density  can  be  carried  out,  using 
Ewald's  technique  for  handling  the "lattice  summations.  This  is  out- 
lined in  Appendix  IV,  the  result  being 

0  y        ^  (2.38) 

Where  Q^  =  J  P(r)   d?  is  that  part  of  the  charge  density  found 
outside  the  spheres,  a=  sV?  R  is  the  lattice  parameter  and  c'=2. 41583 
is  a  constant  characteristic  of  the  fee  structure.  The  superscripts 
on  V^  indicate  the  potential  inside  and  outside  the  spheres. 

Equations  (2.38)  give  us  a  straightforward  way  of  generating 
a  new  potential,  once  the  spherical  average  of  the  charge  density  has 
been  computed  using  the  fonnulas  of , Appendix  III.  This  can  now  be 
used  as  the  input  to  the  energy  band  calculation.  However  this  pro- 
cedure very  often  leads  to  a  divergence,  i.e.,  the  difference  between 
the  potentials  of  successive  iterations  increases.  As  in  atomic  calcu- 
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lations,  it  is  found  that  convergence  is  usually  obtained  by  the 
use  or  an  extrapolation  procedure.  The  simplest  of  these  is  the 
arithmetic  average.  If  (V-jj^,  ^out^  ^^   ^^®  initial  (muffin-tin) 
potential,  and  (Vj^n,  V  +)  is  the  iterated  potential,  then  we  take 
the  new  potential  to  be  (V^^,  V  . )  defined  by 

v,^  =  '/t  [  U '  v..;]  . 

Even  better,  empirically,  is  a  scheme  invented  by  Pratt  -^ ' 
and  used  to  advantage  by  Herman  &  Skillman  -^  in  which  the  results  of 
two  iterations  are  used  to  extrapolate  a  new  starting  potential.  This 
is  best  illustrated  by  a  diagram  (figure  2.2) 


(2.39) 


ITERATED 
''  POTENTl^L 


V    V    V  =  V* 


NITIAL 
POTENTIAL 


Figure  2.2.  Extrapolation  of  a  new  starting  potential 
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In  the  figure,  the  abcissa  represents  the  starting  potential, 
and  the  ordinate,  the  potential  obtained  after  the  iteration  by- 
equation  (2.38).  Each  calculation  is  therefore  represented  by  a  single 
point  on  the  diagram.  The  next  starting  potential  in  this  procedure 
is  chosen  at  the  intersection  of  the  line  joining  the  points  repre- 
senting two  successive  iterations  with  the  45°  line.  If  (V-,  ■->  V-, ) 
and  (V2-^V2)  are  the  initial  and  iterated  potentials  for  the  two 
iterations,  then  the  new  potential  Vo  is  defined  by  the  equations; 

V-   —  rout  *  I I 

Root  -  1 

Pin  -  1 

[V,'-"  -  v;^^]  -  [vj"-^-^j 


(2.40) 

(Since  the  zero-point  of  the  potential  is  arbitrary,  the  extrapola- 
tion procedure  uses  the  potentials  with  respect  to  their  constant 
value  between  the  spheres.) 

It  was  found,  during  this  calculation,  as  Herman  and 

32 
Skill  man  -^  had  for  atoms,  that  quick  convergence  is  obtained  if  the 

two  schemes  (2.39)  and  2.40)  are  alternated. 


CHAPTER  III 
RESULTS  OF  THE  CALCULATION 

3.1.  The  Self -Consistent  Spin-Dependent 
Energy  Bands  in  Nickel; 

The  self-consistent  energy  eigenvalues  were  computed  for 
tvro  values  of  the  exchange  potential:  (i)  the  averaged  free  electron 
value  V^f *  =  -  6(6 p/  Sir)  and  (ii)  the  value  suggested  by  Kohn  and 
Sham   =  2/3  V^^  .  The  initial  potential  chosen  was  a  superimposed- 
atom  potential,  as  defined  in  Section  2.3,  generated  from  self -consis- 
tent atomic  orbitals  corresponding  to  the  configuration 
3^=  5.0,  3d^=  4.0,  4sc(=  4s^=  0.5  . 

The  value  of  the  lattice  parameter  used  in  the  calculation 

38 
is  the  one  quoted  by  Wyckoff,    i.e., 

a  =  6.6586  atomic  units. 
The  radius  of  the  sphere  about  each  lattice  site  corres- 
ponding to  this  value  is  R  =  2,3544  atomic  units.  The  summations  over 
X  in  the  matrix  elements  (AI  -26)  were  taken  up  to  a  value  of  -^ .^^^   6 
and  the  maximum  reciprocal  vector  magnitude  (  vAiich  determines  the 

size  of  the  secular  equation)  was  K    =  6'n"/  a.  By   calculating  the 

max 

energy  eigenvalues  for  a  few  points  for  larger  values,  it  was  found 

that  the  chosen  values,  of  /  ^^  and  K  ^^  were  sufficient  to  insure 
'     max      max 

a  convergence  in  the  energy  values  of  less  than  .005  Ry.  The^ 

"d-like"  states  tend  to  converge  more  slowly  than  the  "s  and  p  -  like" 
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states,  because  of  the  large  negative  values  of  the  logarithmic 
derivatives  for  ^=  2  .  This  implies  that  the  APW  basis  functions 
corresponding  to  the  d-states  vrill  have  a  large  discontinuity  in 
their  derivatives,  and  more  basis  functions  are  necessary  to 
"smooth  out"  the  wave  function. 

The  eigenvalues  and  the  corresponding  wave  functions  were 
calculated  at  the  32  points  of  a  cubic  mesh  in  reciprocal  space 
(grid  #1  of  table  3.1)  during  the  course  of  each  iteration.  Because 
of  symmetry,  actually  only  6  points  were  calculated,  and  the 
corresponding  charge  densities  were  multiplied  with  the  appropriate 
weights  to  derive  the  total  charge  density.  For  the  last  two  itera- 
tions of  each  calculation,  the  mesh  was  reduced  to  half  the  former 
size,  for  an  equivalent  of  256  points  (grid  #2  of  table  3.1). 

It  was  found  that  6  or  7  iterations  were  sufficient  to 
achieve  a  convergence  of  the  energy  levels  to  within  ,01  Rydbergs. 
More  accuracy  than  this  was  deemed  unnecessary,  because  of  the 
uncertainties  in  the  one-electron  potential  and  the  approximations 
inherent  in  the  rauf fin-tin  form  of  this  potential.  The  superimposed- 
atom  potential  used  for  the  first  iteration  turned  out  to  be  quite 
far  from  the  self-consistent  solution,  with  the  result  that  the  first 
two  iterations  appeared  to  be  diverging  (See  figure  3.1).  However, 
^  the  extrapolation  procedure  described  in  Section  2.5  was  sufficient 

to  give  eventually  the  self-consistent  solution. 

In  each  iteration,  after  the  energies  were  calculated,  the 
Fermi  energy  E,  was  estimated  by  counting  the  states  of  both  spins 
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Table  3.1.  Coordinate  grids  in  reciprocal  space. 

Grid  #1  (32  equivalent  points) 
k*(a/ir)     Symmetry^   Weight 


0  0  0 

r 

1 

0  0  1 

A 

6 

0  0  2 

X 

3 

Oil 

S 

12 

0  12 

w 

6 

111 

L 

4 

Grid  #2  (256  equivalent  points) 
ld«-(a/2-rr)    Symmetry   Weight        k*(a/2Tr)    Symmetry   Weight 


0  0  0 

r 

1 

0  2  3 

24 

0  0  1 

A 

6 

0  2  4 

w 

6 

0  0  2 

A 

6 

0  3  3 

K 

6 

0  0  3 

A 

6 

111 

A 

8 

0  0  4 

X 

3 

112 

24 

oil 

2 

12 

113 

24 

0  12 

24 

114 

U 

6 

0  13 

24 

12  2 

24 

0  14 

Z 

12 

12  3 

a 

24 

0  2  2 

2 

12 

2  2  2 

L 

4 

^  The  symmetry  symbols 
Wigner,^° 

are  those  of  Bouckaert, 

Smoluchowski  and 
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in  order  of  increasing  energy  until  10  bands  were  filled,  corres- 
ponding to  the  10  valence  electrons  of  the  nickel  atom.  For  example, 
if  the  energies  had  been  computed  at  the  32  points  of  grid  #1,  then 
the  lowest  320  states  were  chosen  as  the  occupied  levels  for  the 
generation  of  the  next  potential.  Because  of  this  procedure,  the 
configuration  (i.e.,  the  effective  number  of  s,p,  and  d-like  states) 
can  change  from  iteration  to  iteration,  unlike  the  usual  atomic  self- 
consistent  calculation,  in  which  the  configuration  is  held  fixed. 
This  might  be  expected  to  lead  to  a  "collapse"  of  the  energy  bands, 
i.e.,  the  nrmiber  of  occupied  levels  of  the  minority  spin  might  tend 
to  increase  with  each  iteration  until  finally  it  has  been  become 
equal  to  that  of  the  majority  spin.  In  the  course  of  the  calculation, 
it  was  found  that  the  difference  between  the  nimiber  of  electrons  of 
either  spin  did  decrease  slightly,  but  levelled  off  at  a  constant 
value.  This  may  be  a  fortuitous  result  of  using  a  finite  number  of 
points  in  the  Brillouin  zone  to  calculate  the  density  but  the  fact 
that  this  constant  value  is  close  to  the  experimental  is  encouraging. 

The  calculation  was  carried  out  on  an  IBM  709  computer. 
The  time  taken  per  iteration  was  approximately  one  hour  for  grid  #1 
(32  points  in  the  Brillouin  zone)  and  three  hours  for  grid  #2  (256 
points  in  the  Brillouin  zone). 

Figure  3*2  shows  the  final  bands  for  the  first  calculation 
(denoted  by  Vg)  in  •vriiich  the  exchange  potential  was  chosen  to  have 
the  averaged  free  electron  value,   ''^ 
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Here,  the  energy  values  are  plotted  along  the  three  main  symmetry 
directions  in  reciprocal  space: 

A:     r(o,o,o)  -^  x(^,o,o) 

Figure  3 •3  shows  the  same  information  for  the  second  calculation 
(denoted  by  V4),  for  which  V^  =  2/3  V^^^  , 

Both  figures  have  the  same  general  shape,  which  is  essen- 
tially determined  by  the  syrametiy  of  the  face-centered  cubic  lattice. 
Both  show  the  feature,  typical  of  the  transition  metals,  of  a  narrow 
(  '^  .3  Ry.)  "d-band"  imbedded  in  a  much  broader  ('^.Q  I^.)  "s-p"  or 
"conduction"  band.  The  splitting  between  the  two  spin  bands  is  not 
uniform  throughout  the  zone,  but  is  approximately  .07  ^  for  the 
d-like  states  and  zero  in  the  conduction  states. 

The  principal  difference  between  the  two  sets  of  results 
is  in  the  width  and  position  of  the  d-bands.  These  are  considerably 
wider  (-^30^)  sind  higher  ('**.15  Ry.)>  with  respect  to  the  conduction 
band,  for  V^  than  for  Vo  (See  table  3.2).  The  reason  for  this  is  that 

reducing  the  exchcinge  effects  elevates  the  potential  by  an  amount 

1/3 
proportional  to  p  '  .  Therefore,  the  energies  of  the  d-type  electrons, 

vrtiich  are  mostly  confined  to  the  higher-density  interior  of  the  atom, 

are  raised  more  than  the  energies  of  the  conduction  band  electrons, 

irfiich  are  found  further  out  vrtiere  the  density  is  less.  This  elevation 

of  the  d-bands  tends  to  spread  out  the  wave  functions  thus  increasing 
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Table  3. 2. Comparison  of  energy  differences  for  various  potentials. 


Potential 

d-s 

3ep£iration 

d  width 

s-p  width 

splittingi^ 

fli-f-i 

Xy-r, 

Xy-X, 

x;-r, 

.070 

AC, 

Vie<b 

.481^ 

.627 

.320 

.837 

.016 

V/ 

.537 

.687 

.349 

.838 

v,<^ 

.437 

.571 

.299 

.836 

.070 

.013 

V2> 

.491 

.633 

.325 

.837 

¥301^ 

.322 

.441 

.251 

.833 

.076 

-.005 

V3^ 

.403 

.531 

.282 

.834 

V^ 

.478 

.625 

.330 

.841 

.060 

-.004 

V 

.5i^ 

.698 

.362 

.8A2 

V5e(^ 

.488 

.639 

.347 

.865 

.067 

.018 

V5^ 

.535 

.695 

.372 

.866 

The  energy  differences  are  given  in  I^dbergs. 

Vn  =  a  superimposed-atom  potential  with  V^  =  V^^,  corresponding  to 
the  atomic  configuration  (3dc(  =5.0,  3d^  =4.4,  4s 0^=  4sJ  =  0.3). 

°  V2  =  a  superimposed-atom  potential  with  V^  =  V^S  corresponding  to 
the  atomic  configuration  (3do(  =  5.0,  3d^  =  4.0,  4so(=  4s ji  =  0.5). 

V3  =  the  self-consistent  potential  with  Vx  =  V^^ 

®  V^  =  the  self-consistent  potential  with  V^  =  (2/3)  V^*. 

■^  V5  =  Wakoh's  potential  of  reference  10. 
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the  d-band  width. 

Although  the  general  shape  of  the  bands  is  quite  similar, 
the  change  in  the  position  of  the  d-bands  with  respect  to  the  con- 
duction band  is  sufficient  to  alter  some  of  the  details.  For  example, 
at  the  L  -  point  (  TT  /a,  Ti'/a,  Tf/a)  there  are  two  doubly  degenerate 
d-like  levels  Lj,  and  two  singly  degenerate  levels:  In  (  a  mixture  of 
s  and  d-like  states  )  and  I^  (which  has  p-like  symmetry).  Figure  3.4 
shows  the  bands  in  the  region  around  this  point . A  is  the  line  from 
r (0,0,0)  to  L,  and  Q  is  the  line  perpendicular  to  A  on  the  face  of 
the  zone,  going  from  L  to  W  (Tr/a,2tr/a,  0  ).  It  can  be  seen  from  this 
figure  that  the  energy  bands  are  substantially  different  in  the  two 
cases,  the  p-like  level  Lj  being  much  lower  with  respect  to  the 
d-bands  for  Mj^  than  for  V^.  This  is  particularly  significant  in  that 
the  position  of  the  Fermi  level  with  respect  to  the  conduction  band, 
and  consequently  the  Fermi  surfaces,  are  considerably  altered  by  the 
shifting  of  the  I^  -  level.  Table  3.3  lists  the  ordering  and  the 
resultant  Fermi  surfaces,  and  shows  that  the  reduction  of  the  exchange^ 
effects  has  qualitatively  changed  these  surfaces. 

The  experimental  data  (See  Section  3.3)  favor  the  second 
calculation  (V^),  because  of  the  multiply-connected  Fermi  surface 
in  the  sijcth  c<-band.  This  surface  is  simlar  to  the  famous  copper 
Fermi  surface,  having  "necks"  along  the  (1,1,1)  direction. 

3.2.  Comparison  with  Previous  CalcTilations ; 
Within  the  last  five  years,  several  calculations  or  estimates 
of  the  energy  band  structure  of  nickel  have  appeared  in  the  literature, 
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Figure  3.4.  Detail  of  the  energy  bands  at  an  L-point 
for  the  two  self-consiBtent  potentials 


Table  3.3.  Detail  at  an  ]>point. 


Potential       Ordering  at       Resultant  shape  of  Fermi 
the  L-point       surfaces  near  the  L-point 


V3  <^  L3  <  Ef  <  Lj        Closed  surface  in  the 

6th  band. 


^3  -/*  -^f  "^  ^  "^  ^        Closed  surfaces  in  the 

5th  and  6th  bands,  holes 
in  the  4th  band. 

\°^  L2<L2<Ef        Multiply-connected  surf- 

ace in  the  6th  band  with 
a  "neck"  at  the  L-points, 

^U-^  'h.^^i'^^  Closed  surface  in  the 

6th  band. 
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The  augmented  plane  wave  method  has  been  applied  to  the  paramagnetic 
case  by  (i)  Hanus,   who  used  a  potential  generated  from  renormalized 
atomic  orbitals  of  the  configuration  (3d)^(4s)^,  (ii)  Mattheiss,^ 
Tidio  used  a  super imposed-atom  potential  corresponding  to  the  atomic 
configuration  (3d)9(4s)^,  (iii)  Snow,  Wftber  and  Switendick,  ^  who 
tried  several  different  superingaosed-atom  potentials  with  configu- 
rations (3d)-'-°"^(4s)^,  0  ^  X  ^  2. 

The  three  calciilations  are  qualitatively  similar  to  each 
other,  differing  only  in  the  position  of  the  d-bands  with  respect 
to  the  conduction  band.  In  reference  7  it  is  shown,  as  confirmed  by 
the  results  presented  in  Section  3.1,  that  this  d-band  shift  has  a 
profound  effect  on  the  topology  of  the  Fermi  surface,  caused  mainly 

by  the  altered  ordering  of  the  L  -  point  levels.  (See  table  3.4.) 

8""10 
Yamashita  and  Wakoh,     have  also  done  several  calculations 

for  nickel,  both  in  the  paramagnetic  and  ferromagnetic  states,  using 

the  Green's  function  method.  Since  this  method  is  formally  equivalent 

to  the  AFW  method,    their  results  are  slLso  very  similar,  showing 

the  same  sensitivity  of  the  d-bands  to  changes  in  the  potential. 

In  order  to  compare  all  of  these  results,  we  note  that  the 

results  of  each  calculation  can  be  fairly  well  described  by  only  two 

parameters : 

(i)  the  position  of  the  d-band  with  respect  to  the  conduction 

band,  which  can  be  characterized  by  the  energy  difference  between  the 

states  125  ^^^    '  t  and 

(ii)  the  d-band  width,  a  measure  of  ii^ich  is  the  energy 

difference  between  Xc  and  1-, . 
p     1 
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Table  3.4.  Results  of  Snow,  Waber  and  Svdtendick? 
for  paramagnetic  nickel 


Atomic  configuration  Ordering  of 
used  for  the  super-  energies  at 
imposed-atom  potential   the  L-point 


(3d)8.0  (4s)2.0      L3^E^<I^ 


(3dr-5  (^3)1.5      I3<E^<I^ 


Fermi  surface  in 
the  6th  band 


Closed  surface  with  pro- 
trusions toward  the  X- 
points. 


(3d)9-0  (^^ji.o 


(3d)9-5  (^gjO.5 


I^<L3 


<E. 


I^<Ejf<l3 


(3d)^°'°  i<^<Er<L3 


Multiply-connected  surf- 
ace, with  "necks"  at  the 
L-points . 

Closed  surface  with  pro- 
trusions toward  the  X 
and  L-points. 


L 
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Using  these  two  parameters  as  coordinates,  each  calculation 
is  plotted  on  the  graph  in  figure  3.5.  For  comparison,  a  few  calcu- 
lations for  other  elements  (Cu,Co,Fe)  with  an  fee  structvire  are 
included.  We  note  that  no  matter  what  the  potential,  ^^"11  the  calcu- 
lations tend  to  lie  on  a  straight  line,  which  suggest  that  the  ener- 
gy band  structure  of  an  fee  transition  metal  could  be  described 
by  a  single  parameter,  say,  the  position  of  the  d-bands. 

This  diagram  shows  clearly  the  considerable  variation  between 
the  resxilts  of  previous  work.  The  principal  differences  can  be  traced 
back  to  the  atomic  configuration  used.  Those  calculations  which 
assumed  a  smaller  number  of  d-eleetrons  are  fovmd  toward  the  lower 
left  hand  comer  of  the  diagram.  Although  all  of  the  calculations 
listed  here  used  the  averaged  free-electron  exchange  potential  V^  , 
the  effect  of  reducing  the  exchange  moves  the  results  in  the  opposite 
direction. (See  Section  3.1.)  In  short,  anything  which  makes  the 
potential  more  attractive  for  the  d-eleotrons  narrows  the  d-band 
and  lowers  it  with  respect  to  the  s-p  band. 

Of  the  previous  calculations  on  transition  elements,  only 
two  have  been  cai'i'ied  to  self-consistency,  viz.,  that  of  Wakoh 
on  copper  and  ferromagnetic  nickel,  and  that  of  Snow  and  Waber 
on  copper,  both  of  which  used  V^  .  We  can  see  by  looking  at  figure 
3.5,  that  there  is  a  considerable  difference  between  the  tvro  results 
for  Cu  (references  10  and  41)  although  both  lie  near  the  straight 
line,  Wakoh 's  d-bands  are   displaced  upward  from  those  of  Snow  and 
Waber  by  approximately  .15  Ry.  This  is  almost  exactly  the  same  as  the 
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•      ITERATION  OF    THIS   CALCULATION 
D     PREVIOUS    CALCULATIONS 


INCREASINO 
NUMBER  OF 
D  ELECTRONSi 


INCREASING 

ATOMIC 

NUMBE 


NI-MATTHEISS(6) 


CU-MATTHE1SS(6)D 
CU-ARLINGHAUS  \(41D 

rcu-BURDicK  -Vr^ 

(16)  '-'° 

CU-WAKOH(IO) 


0.3 


0.2 


Nip- WAKOH  (10) 

C0-MATTHEISS(6) 
D 

NIB-V  FE(*CC)-W00D  (14) 

NIei-WAKOH  (10) 
Nla-V4 

D— NI-HANUS  (5) 

Nl-  YAMASHITA    a    WAKOH  (9) 

NI-SNOW     ET  ALC?) 
Nl^-V3 


CU-SEGALL  (15) 
Nla-V3 


CU-SNOW   a  WABER    (41) 


0.2  0.3  0.4  0.5 

D-BAND  WIDTH  =  E(X5)-E{X,)RYDBER0S 

Figure  3.5.  Comparison  of  transition  element 
energy  band  calculations 
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displacement  between  Wakoh's  bands  for  ferromagnetic  Ni  and  those 

—■fe 
of  V^  (the  self -consistent  bands  for  V^^  ).  The  reason  for  this 

discrepancy  lies  in  the  simplified  version  of  self-consistency  used 
by  Wakoh.  The  wave  functions  used  to  generate  a  new  potential  after 
each  iteration  were  not  chosen  uniformly  over  the  BriUouin  zone.  His 
procedure  was  to  pick  out  5   fionctions  representative  of  the  d-like 
electrons  (corresponding  to  local  maxima  in  the  density  of  states 
curve)  and  another  function  representing  a  conduction  electron, 
corresponding  to  energy  E(f^)  +  (3/5)E^.  These  fimctions  were  then 
weighted  to  give  5.0  do^-electrons,  4.4  d_p -electrons  and  0.3  conduc- 
tion electrons  of  both  spins.  For  the  copper  calculation,  the  corres- 
ponding weights  were  10  d-electrons  and  one  conduction  electron.  Since 
the  "conduction  electron"  function  chosen  here  is  a  mixture  of  s  and 
d-like  ftmctions,  this  procedure  tends  to  fix  the  amount  of  d-like 
character  at  too  high  value,  thus  displacing  the  d-bands  upward.  His 
final  result  is  very  close  to  that  of  V,  (the  self -consistent  bands 
for  2/3  V^c  )»   since  the  effect  of  overestimating  the  exchange  tends 
to  balance  the  effect  of  using  too  many  d-electrons. 

For  comparison,  the  effective  number  of  d-lilce  electrons 
calculated  for  Ni  from  the  self -consistent  results,  according  to 
the  formulas  of  Appendix  III,  are  listed  in  table  3.5.  We  note  that, 
although  the  effective  number  of  d-electrons  is  about  the  same  (-»'  9,4) 
for  the  superimposed-atom  potential  V  and  the  self-consistent  poten- 
tial V3,  the  resulting  bands  are  displaced  '^,11^,  from  each  other. 
This  is  a  reflection  of  the  distortion  of  the  d-functions  on  going 
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Table  3.5.  Charge  distributions. 


Potenticil 

Outside 
charge 

Inside  charge 

Total 
charge 

Effective 
d-charge 

Qo 

Qs 

Qp 

Qd 

\ 

0^^ 

V30I 

.26 

.20 

.14 

4.72 

5.33 

4.97 

V3J3 

.28 

.20 

.16 

4.02 

4.67 

4.38 

V40t 

.35 

.24 

.24 

4.47 

5.31 

4.79 

\fi 

.39 

.24 

.26 

3.79 

4.69 

4.14 

Cu-SAP^ 

.58 

.50 

.AA 

9.45 

10.97 

9.96 

Cu-B 

.61 

.51 

.50 

9.33 

10.95 

9.90 

QqI  =  Qd(Qt/(Qt~Qo))'  The  quantities  Qq  etc.  are  defined  in  Appen- 
dix III. 

The  results  for  Cu  are  taken  from  Arlinghaus^  vrfio  calculated  the 
bands  for  a  superimposed-atom  potential  (Cu-SAP)  and  for  the  poten- 
tial used  by  Burdickl°  (Cu-B). 
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to  self-consistency, The  atomic  d-functions  are  more  localized  than 
the  self-consistent  ones,  which  has  the  effect  of  elevating  the 
bands  through  the  Cotilomb  repvdsion  term.  There  is  a  similar  displa- 
cement between  the  self-consistent  copper  bands  of  Snow  and  Waber  ^ 
and  those  of  Burdick    and  Arlinghaus,  ^  even  though  the  nranber  of 
d-electrons  is  close  to  10  in  both  cases.  (See  table  3.5.) 

3*3.  Comparison  With  Experimental  Data 
Of  the  wealth  of  experimental  data  that  has  been  published 
on  nickel,  there  are  several  firmly  established  facts  which  should 
be  e3q)lained  by  an  energy  band  calculation,  if  the  model  is  to  have 
any  validity  at  all:  (1)  the  saturation  value  of  the  magnetization, 
(2)  the  electronic  specific  heat  at  low  temperatxa'es,  (3)  the  Fermi 
surface  topology  as  deduced  from  de  Haas-van  Alphen  measurements, 
(4)  the  saturation  of  the  magnetoresistance,  and  (5)  the  negative 
spin  density  found  from  neutron  diffraction  and  positron  annihilation 
data. 

In  order  to  interpret  the  first  two  of  these  properties,  we 
need  to  calculate  from  the  energy  bands,  the  density  of  states,  i.e., 
the  number  of  allowed  energy  levels  per  unit  energy.  Mathematically 
this  is  '''''■ 


(T.yisJ.^      ' 


-» 


vAiere  Eg(  k  )  is  the  energy  of  an  electron  with  spin  s  and  reduced 

wave  vector  k  and  12  is  the  volume  of  the  unit  cell.  Since  E  is  not 

s 
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known  analyticalily  from  the  computation,  we  approximate  (3.1)  hy 
.^Ce)    -    ^^     Z    ^/^^  (3.2) 

where  the  sum  is  over  those  t  points  at  vrfiich  Eg  ( "k  )  is  known,  and 
w  (  k  )  is  the  weight  for  the  point  ^.»  (See  table  3.1.) 

In  this  calculation,  E^  (  if )  was  calculated  at  only  20 
non-equivalent  points  of  the  Brillouin  zone,  i^ich  is  not  enough  to 
give  an  accurate  n^  (E).  In  Appendix  V  is  shown  a  procedure  for  inter- 
polating Eg  (  k  )  at  enough  points  to  give  sufficient  accuracy. 

The  calculated  ng  (E)  for  both  spins  is  shown  in  figure  3.6 
for  the  two  self-consistent  potentials  V3  and  V^.  It  can  be  seen  that 
the  narrowing  of  the  d-bands  in  V3  substantially  increases  the  height 
of  the  peaks  in  ng  (E). 

The  Fenni  level  E^  is  determined  by  the  relation, 

vrfiere  n^  is  the  number  of  valence  electrons,  10  in  nickel,  and  I^ 
is  the  minimum  energy  of  the  corresponding  bands.  The  difference 
in  the  number  of  electrons  of  either  spin  is  then 


r  ^ 


(3.4) 

•min 

The  measured  value  of  yi*  corresponding  to  the  saturation 
value  of  the  magnetization  is  O.6O6  electrons/atom  (reference  43 
page  317)  . 
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Energy  (Rydbergs)  with  respect  to  VqO( 
Figure  3.6.  The  density  of  states  curves 
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The  calculated  values  ofyu.  for  the  super imposed-atom 
potentials  using  V^  are  slightly  larger  than  the  esqperimental 
value,  ■which  tends  to  substantiate  the  conjecture  that  this  appro- 
ximation overestimates  the  exchange  effects.  Going  to  self-consis- 
tenty  slightly  reduces  the  value  bringing  it  closer  to  ,6,  i.e., 
/(.(Vo)  =  .65,  a(V,)  =  .62,  The  error  in  the  calculation  of  yn  is 
fairly  large, A^-  np(E^)AE^,  since  n^  (E^)  is  high  (~40  electrons/ 
atonw-I^.)  An  error  of  AE^  "J, 001  Ry.  therefore  gives  an  error  in  ^ 
of  "^ .04  electrons/atom,  so  that  the  agreement  here  is  quite  satisfac- 
tory. 

The  electronic  specific  heat,  c  ,  at  low  temperatures  is 
related  to  the  density  of  states  at  the  Fermi  energy  by  the  equation 


c,  =    ^*  [\C^/)  ^  ^j(^^^] 


(3.5) 
3 

(reference  45,  page  150).  The  measured  value  of  c^  corresponds  to  a 

total  density  of  3.1  electrons/atom-ev.  The  calculated  values  of 

n(E£)  are  n(Vo)  ^   4.5j  n(V,  ) -^  3.0  electrons/atom-ev.  The  discrepancy 

in  the  value  for  Vo  is  a  further  indictment  of  these  bands,  since  the 

method  of  calculating  n(E)  should  be  expected  to  give  a  density  lower 

than  the  exact  result. 

Both  calculations  show  the  Fermi  energy  occurring  at  a  high 

peak  of  the  density  of  states  curve  Tirfiich  agrees  with  Slater's 

criterion,  ^'  for  ferromagnetism.  The  results  are  also  confirmed  by 

the  recent  measurements  of  the  magnetic  susceptibility,  ^  X,  ,  irfiich 

satisfies  a  relationship  of  the  form; 
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The  measurements  indicate  a  very  small  Xj,  vrtiich  would  be 
obtained  only  if  n^^  (E^ )  is  small.  We  note  that  in  both  Vo  and  V, , 
since  the<^d-bands  are  filled,  n^^  (E^)«ntt(E^). 

Although  the  Fermi  surface  of  nickel  is  not  as  completely 

mapped  out  as  that  of  Cu  and  the  noble  metals,  at  least  one  featiire 

k9   50 
has  been  firmly  established.  The  de  Haas-van  Alphen  measurements,  ' 

indicate  a  small  cross-sectional  area  of  the  Fermi  surface  in  the 
(1,1,1)  direction.  This  is  interpreted  as  a  "neck"  in  a  Fermi  surface 
similar  to  that  found  in  copper,  i.e,,  roughly  spherical  with  protru- 
sions that  make  contact  with  the  Brillouin  zone  at  the  L-points,  The 
neck  area  is  about  l/lO  that  found  in  Cu,  corresponding  to  an  angle 
of  6.8°  +  0.2°  subtended  at  the  P  -point.  ^°  If  the  energy  surfaces 
are  assumed  to  be  hyperbolic, 

£  =  JL 1    -  _-?  (3^7) 

then  m^  and  m^  are  the  transverse  and  longitudinal  effective  masses, 


v^ose  measured  values  are 
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m^  =  ,26  +  .04,  flijj  =  ,65  +  .10 
The  Fermi  surface  derived  from  the  energy  bands  for  Vi  is 
shown  in  two  cross-sections  in  figure  3.7.  It  consists  in  this  case 
of  four  sheets  for  the  down-spin  (fl )  bands  and  one  in  the  up-spin 
( o( )  bands .  The  ^  -   sheets  consist  of  small  hole  surfaces  at  the 


58 


> 

u 
o 


•p 

o 

n 
I 

a 

ra 
o 

u 


59 


X-points  in  the  3rd  and  4th  bands,  a  closed  surface  with  protrusions 
along  the  Z.   directions  in  the  5t.h  band,  and  closed  surface  with 
protinisions  along  the  A  -directions  in  the  6th  band.  The  sole 
o(  sheet  is  a  multiply-connected  surface  with  necks  at  the  L-points, 
in  qualitative  agreement  with  experiment.  The  calculated  angle  sub- 
tended at  r  is  10°,  slightly  larger  than  the  measured  value.  The 
theoretical  values  of  the  effective  masses  at  the  neck  are  in  much 
worse  agreement;  m^  =  .5,  m_g  =  1.1. 

This  discrepancy  is  probably  due  to  the  fact  that  the 
energy  stirfaces  at  lyj  do  not  conform  to  the  simple  equation  (3.7)> 
but  is  more  like  the  solution  of  a  quadratic  equation 

'A 

£=  -^(s^cLk')    ^  i  [(Si^ak')"-    ^(/>^/>k')]         (3.8) 

along  a  particuliir  direction  in  K-space,  vrtiere  s  =  ECLo)  +  E(L2)  and 
p  =  E(L3)-  E(I^). 

The  principal  result  of  this  behavior  is  that  the  effective 
mass  m   =(3E/  ^k)^,-   varies  strongly  with  the  position  of 
£_£■,  unlike  (3.7)  where  it  is  a  constant.  In  fact,  shifting  the  level 
Lp  up  with  respect  to  E^  by  an  amount  of  '«  .02  Ry  would  bring  the 
effective  masses  into  approximate  agreement  with  the  experimental 
values.  A  shift  of  this  magnitude  would  also  bring  the  neck  size  into 
agreement  with  the  measured  value. 

The  model  of  the  ol   Fermi  surface  is  consistent  with  the 
magnetoresistance  date  of  Fawcett  and  Reed.  ^^»^^   These  measurements 
clIso  show  that  the  magnetoresistance  saturates  at  high  magnetic  fields. 
This  type  of  behavior  is  usually  associated  with  an  uncompensated 
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material.  Nickel  being  an  even-valenced  metal,  would  ordinarily  be 
e^qjected  to  behave  like  a  compensated  material.  ^^  This  would  be  the 
case  in  the  ordinary  energy  band  model  in  vdiich  the  two  spins  are 
degenerate.  However,  the  unrestricted  case  can  give  rise  to  the 
situation  where  a  different  number  of  bands  are  occupied  for  either 

spin,  the  sum  of  the  occupation  numbers  being  an  odd  number. 

51 
ExperimentaUy,  ''     the  number  of  electrons  per  atom  must  satisfy  the 

relation 

-y^^U)    ■*■   -^eCf)  -   ^f-^^-^yj^^J  =10   (3.9) 

where  n^  (s)  is  the  number  of  elecbrons/atom  on  the  electron  Fermi 
surfaces,  and  nj^  (s)  is  that  on  the  hole  Fermi  surfaces.  Equation 
(3.9)  is  satisfied  by  both  of  the  energy  bands  Vo  and  V,,  since  in 
both  cases  five  of  the  o{  -bands  and  four  of  the  J3  -bands  are  filled, 
so  that  dnOy  one  of  the  ten  valence  electrons  is  left  to  be  distri- 
buted over  the  Fermi  surfaces. 

Recently,  it  has  been  discovered  through  measurements  of 
positron  annihilation  ^-^  and  neutron  diffraction  ^^  in  ferromagnetic 
Mi,  that  the  electron  spin  density  is  negative  in  the  outer  regions 
of  the  unit  cell.  It  is  found  that  for  both  V3  and  V, ,  the  conduction 
electrons  are  polarized  opposite  to  the  d-electrons,  so  there  is  a 
theoretical  negative  spin-density  outside  the  spheres,  and  just  inside 
the  spheres.  The  measured  ^4  value  of  19^  of  the  magneton  number  /*- 
(.606)  is  larger  than  the  theoretical  ones  for  both  calculations, 
perhaps  because  of  the  neglect  of  the  non-spherical  terms  of  the 
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potential.  As  pointed  out  in  reference  4,  these  non-spherical  terms 
are  essential  for  the  interpretation  of  the  neutron  diffraction  data, 
and  will  have  to  be  included  before  quantitative  agreement  can  be 
obtained. 

There  are  also  a  great  deal  of  optical  data  available  for 
nickel.  There  are  summarized  in  references  1  and  2,  which  analyze 
the  ferromagnetic  Kerr  effect,  the  dielectric  constant  and  the 
conductivity.  The  structure  found  in  these  data  (at  .3,  .6  and  1.4  ev) 
is  interpreted  as  interband  transitions.  Because  of  the  errors  involved 
in  calculating  a  density  of  states  curve  (See  Appendix  V),  it  is 
difficult  to  determine  i^ether  or  not  structure  on  this  small  a 
scale  is  present.  It  can  be  seen  from  figures  3.2  and  3.3  that 
transitions  to  empty  states  above  the  Fermi  level  are  possible  in  the 
regions  around  the  X  and  L  points.  These  are  the  most  probable 
since  the  density  of  states  is  high  at  these  points.  However,  the 
density  of  states  curve  calculated  here  is  not  accurate  enough  to 
compare  the  theoretical  energy  bands  with  the  experimental  optical 
data. 

A  summary  of  the  ejqierimental  and  calculated  data  is  found 
in  table  3.6.  It  can  be  seen  that  the  self-consistent  bands  for 

^x  ^  ^/^  ^x  ^^4^  ^®  ^  ^'^^^   better  agreement  than  those  for 
V  =v^(Vo) 
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Table  3.6,  Comparison  of  calculated  and  eaqjerimental  quantities. 


Ceilculated      Experimental 
V3       V4 


Saturation  moment 
(Bohr  magnetons) 

Density  of  states  at 
the  Fermi  level 
(electrons/atom-ev) 

Neck  of  the  6<*  Fermi 
surface  (angle  subtend- 
ed at  the  r  -point) 

Effective  masses  at  the 
neck  (in  units  of  the 
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CHAPTER  IV 
CONCLUSIONS 

The  purpose  of  this  work  has  been  to  solve  as  accurately 
as  possible  within  the  limits  of  present  computational  techniques, 
the  unrestricted  Hartree-Fock  (UHF)  equations  in  a  ferromagnetic 
solid.  Whether  or  not  this  is  a  valid  model  will  not  be  discussed  here. 
To  really  test  this  method  would  involve  a  determination  of  the  excited 
state  energies,  which  would  permit  an  analysis  of  the  temperature- 
dependent  properties  in  a  way  similar  to  that  of  the  Heisenberg 
model.  However,  this  calciilation  has  been  limited  to  the  ground  state 
of  ferromagnetic  nickel.  This  follows  naturally  from  the  use  of  the 
energy-band  model,  irfiich  although  it  can  explain  certain  temperature- 
independent  properties,  has  not  yet  been  extended  in  a  rigourous  way. 

In  attempting  to  solve  the  UHF  equations,  we  have  used  the 
one-electron  approximation.  This  has  had  considerable  success  for 
atomic  systems,  and  it  can  be  shown,  ^°  that  it  can  give  a  solution 
that  is  very  nearly  as  good  as  the  exact  Hartree-Fock  value.  There 
is  no  reason  to  believe  that  the  same  should  not  be  true  for  a 
solid  .  The  problem  here  is  to  find  the  "best"  one-electron  approxi- 
mation, i.e.,  to  find  an  effective  local  exchange  potential  that  can 
accurately  reproduce  the  exact  Hartree-Fock  solution.  One  way  to  do 
this  would  be  to  calculate  the  totsil  energy  of  the  solid  in  the 
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correct  vray,  viz.,  E^  =  <'^\ii\'^},  vhere   t|/  is  a  determinant  made  up 
of  the  approximate  self-consistent  one-electron  fvmctions  and  ^ 
is  the  exact  Hamiltonian^  and  then  to  minimize  E,  with  respect  to 
the  exchange  approximation.  lindgren  ^  has  found  that  this  method 
works  very  well  for  atoms,  giving  feetter  agreement  with  the  experi- 
mental binding  energies  than  the  exact  Hartree-Fock  method.  Such  a 
scheme  is  of  course  much  more  difficult  in  a  solid,  but  de  Cicco  ■'•''' 
has  recently-  shown  that  it  is  feasible  to  compute  the  total  energy 
of  a  solid,  so  that  an  investigation  of  this  type  might  be  possible 
in  the  future. 

This  has  not  been  done  in  this  work,  but  is  has  been  shown 
that  slightly  changing  the  one-electron  Hamiltonian  critically 
changes  the  resulting  energy  band  structure.  In  particular,  it  turns 
out  that  using  the  averaged  free-electron  exchange  approximation 
Vx  -  -6(6^ /STr)    gives  qualitative  disagreement  with  experiment, 
vdiereas  reducing  the  exchange  to  (2/3 )Vj  gives  agreement  for  most  of 
the  experimental  data.  The  strongest  argument  for  this  conclusion  is 
the  Fermi  surface  structure  deduced  from  the  calculated  energy  bands. 
The  existence  of  a  copper-like  Fermi  surface  with  a  small  "neck"  in 
the  (1,1,1)  direction,  which  has  been  firmly  established  ^^-52  ^^ 
de  Haas-van  Alphen  and  magneto-resistance  experiments,  is'  not  predict- 
ed  by  the  V^  -bands,  but  does  occur  itrtien  the  exchange  is  reduced. 
This  tends  to  confirm  the  conclusion  (See  Section  2.2.)  that  others 
have  reached  for  atomic  systems,  i.e.,  that  V*^over estimates  the 
exchange  effects. 
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This  sensitivity  of  the  energy  bands  to  changes  in  the 
•one-electron  approximation  was  not  pointed  out  in  previous  calculations 
of  this  type,  mainly  because  they  were  not  carried  out  to  self- 
consistency,  and  therefore  no  conclusions  could  be  made  about  the 
adequacy  of  the  Hamiltonian.  There  has  been  one  other  calculation 
(that  of  Snow  and  Waber  ^  on  Cu)  that  was  self-consistent  using  the 
V  exchange.  The  discrepancy  of  the  Fermi  surface  did  not  occur  in 
this  case,  since  the  position  of  the  Fermi  level  with  respect  to  the 
position  of  the  d-bands  is  different  for  Cu  than  for  Ni.  In  Cu,  there 
are  11  valence  electrons  per  atom,  10  of  which  are  "d-like",  so  that 
the  d-bands  are  full.  The  Fermi  level  therefore  lies  above  the 
d-bands  in  Cu,  whereas  in  Ni,  having  one  less  valence  electron,  has 
its  Fermi  level  slightly  below  the  top  of  the  d-bands.  In  both  cases 
(this  work  and  reference  a),  going  to  self-consistency  lowers  the 
d-bands  with  respect  to  the  s-p  bands.  For  Cu,  the  position  of  E^ 
with  respect  to  the  s-p  bands  remains  essentially  the  same,  but  for 
Ni,  E„  is  lowered  with  respect  to  the  s-p  bands.  As  shown  in  Section 
3.1,  the  position  of  E^  with  respect  to  the  s-p  bands  is  \ihat 
determines  the  shape  of  the  Fermi  surfaces.  Therefore,  going  to  self- 
A        consistency  in  Cu  showed  no  basic  change  in  the  Fermi  surface,  as  is 
the  case  for  Ni.  This  sensitivity  makes  Ni  a  better  test  case  than 
Cu  for  testing  the  validity  of  an  exchange  approximation  such  as  v^  . 

Actually,  before  we  can  definitely  say  whether  the  discre- 
pancy is  due  to  the  inaccuracy  of  V^*,  we  should  examine  the  approxi- 
mations made  in  this  computation,  in  order  to  see  if  there  could 
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possibly  be  some  other  effect  \rfiich  would  shift  the  bands  enough  to 
bring  the  V^-  bands  into  agreement  with  experiment.  From  Figiu-e  3.2, 
it  can  be  seen  that  this  shift  wo\ild  have  to  be  approxunately  .1  Ry, 
in  order  to  bring  the  d-bands  up  far  enough  to  give  the  experimentally 
observed  Fermi  surface. 

There  are  four  important  approximations  which  have  been 
used  in  this  computation: 

(1)  the  neglect  of  relativistic  ef feet s .  They  can  be  handled 
within  the  framework  of  the  energy  band  model,  and  are  certainly 
important  for  heavy  atoms.  Calculations  which  have  been  made  (e.g., 
reference  36)  show  that  these  effects  are  on  the  same  order  of  magni- 
tude in  a  solid  as  in  a  free  atom.  For  the  Ni  atom,  the  relativistic 
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splittings  of  the  one-electron  energy  values  are  -^    ^  .01  Ry.,  so  that 

it  does  not  appear  that  they  could  qualitatively  change  the  Fermi 
surface  structtire. 

(2)  the  neglect  of  the  non-spherical  component  of  the 
potential.  In  deriving  a  potential  from  a  muffin-tin  charge  density 
(See  Appendix  IV)  we  used  the  spherical  average  instead  of  the  actxial 
Ewald  potential.  In  Reference  18,  it  is  shown  how  these  effects  can 
be  taken  into  account,  by  expanding  the  non-muffin-tin  part  of  the 
potentisil  in  terms  of  a  Fourier  series.  For  example,  it  turns  out 
that  the  lyj  level  can  be  well  approximated  by  a  single  symmetrized 
plane  wave  corresponding  to  k  =  (ir/a)  (  1,1,1).  The  Fourier  component 
of  the  Ewald  potential  corresponding  to  k  =  (2Tr/a)  (1,1,1)  is 
equal  to  -(o/a)  (.0096).  Substituting  in  the  values  of  Q  =  .26  and 
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a  =  6.65  shows  that  the  energy  shift  of  the  I^  -  level  would  be 
on  the  order  of  only  .001  Ry.  Since  this  level  is  the  one  which 
determines  the  neck  of  the  Fermi  svirface,  it  appears  that  the 
inclusion  of  these  non-muffin-tin  terms  would  change  the  resultant 
Fermi  surface  only  by  a  negligible  amount. 

(3)  the  neglect  of  the  non-spherical  components  of  the 
charge  density.  The  Ewald  potential  was  derived  on  the  assumption 
of  a  muffin-tin  charge  density.  However,  the  departxires  of  the 
actual  density  from  this  approxijnation  are  certainly  significant, 
as  can  be  seen  from  the  neutron  diffraction  data.  ^  These  effects 
are  much  harder  to  estimate,  ■^'   especially  the  non-spherical  parts 
within  the  spheres.  It  can  be  done  by  an  expansion  in  terms  of 
spherical  harmonics,  but  if  more  than  a  few  -^-values  are  important, 
the  number  of  integrals  to  be  calculated  becomes  prohibitively  large. 
They  should,  of  course,  be  investigated  if  the  data  are  to  be 
accurately  interpreted,  but  it  is  unlikely  that  they  will  seriously 
affect  the  energy  band  structtore  as  presented  here. 

(4)  the  neglect  of  the  distortion  of  the  core  states.  In 
this  calculation,  the  core  state  wave  function  (ls,2s,2p,3s  and  3p) 
were  assumed  to  remain  unchanged  from  the  values  fovind  for  the  Ni 
atom.  Snow  and  Waber  ^  in  their  calculation  on  Cu  tried  to  estimate 
the  error  involved  in  this  asaimption.  They  used  two  Vcilues  of  the 
core  density,  one  obtained  from  an  ordinary  self -consistent  atomic 
calculation,  and  the  other,  from  an  atomic  calciolation  in  which  all 
of  the  wave  functions  were  constrained  to  be  in  the  Wigner-Seitz 
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sphere  (i.e.,  that  sphere  with  volume  equal  to  a  unit  cell).  It  was 
found  that  constraining  the  core  states  had  the  effect  of  pushing 
the  d-states  u£  with  respect  to  the  s-p  band  by  approximately  .03^y. 
During  the  course  of  this  work,  the  bands  corresponding  to  the  3s 
and  3p  levels  were  computed  using  the  actual  crystalline  potential, 
and  a  new  potential  derived  from  the  resultant  wave  functions.  The 
changed  potential  shifted  the  d-bands  up  by  an  amount  less  than  .01 
Ry.  Neither  estimate  is  enough  to  qualitatively  change  the  Fermi 
surfaces,  although  a  more  accurate  calculation  should  certainly 
take  the  distortion  of  the  core  states  into  account. 

The  arguments  presented  here,  although  not  conclusive, 
tend  to  suggest  that  any  discrepancy  found  can  be  attributed  to  the 
approximation  made  in  the  Hamiltonian,  rather  than  in  the  approxima- 
tion made  to  find  the  eigenvalues  of  this  Hamiltonian. 

To  the  list  of  approximations,  coiold  possibly  be  added 
another,  i.e.,  the  neglect  of  any  correlation  effects.  These  are  of 
course  not  included  in  the  Hartree-Fock  model,  but  they  may  be  im- 
portant in  ferromagnetic  nickel.  However,  the  importance  of  correla- 
tion cannot  be  determined  if  the  exchange  is  not  known  accurately. 
Indeed,  it  may  be  that  the  averaged  free-electron  approximation  is 
a  better  estimate  of  the  Hartree-Fock  exchange  effects  than  it  appears 
here,  and  that  reducing  the  magnitude  of  the  exchange  merely  introdu- 
ces some  effective  correlation.  More  study  is  needed  on  these  effects, 
if  the  energy  band  model  is  to  be  successful  in  the  explanation  of 
ferromagnetic  effects. 
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In  conclusion,  the  main  result  of  this  work  is  that  the 
unrestricted  Hartree-Fock  equations  can  be  solved  in  a  solid,  at 
least  in  the  one-electron  model,  and  that  this  scheme  forms  a 
reasonable  model  for  a  ferromagnetic  solid.  The  accuracy  obtained 
is  sufficient  to  give  at  least  qualitative  agreement  with  experiment, 
and  could  probably  be  improved  if  better  approximations  to  the 
Hamiltonian  were  used. 


APPENDIX  I 

*  APW  MTRIX  EUJMENTS 

A  single  APW  basis  function  has  the  form,  for  a  crystal  with 
one  atom  per  unit  cell 

^        /  «»/  ^       _»  (Al-l) 

^  £C/t->-;  4  Z,  "^JL^^^^  '^J'^.^^  yjiJ^^h 

irfiere      E  (k)         is  a  step  function, 

£  and         ('-O^)''    ^      iXL^i)    i.    V(^))(f-kJ=   £(fU^)AAl'3) 

The  spherical  harmonics  ^j^are  chosen  according  to  the  phase 
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convention  of  Condon  and  Shortley,    i.e.. 


r 


so  that  they  have  the  follovring  properties: 


f^L^f^)      'fx'r.'C^)    '^^^      '     ^M'^"'^      . 


(Al-5) 


(Al-6) 


-«    ;r 
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<v.  JL 


e  ""■  -  -     ^^  2,X'  ^^^  ^'^^  *'^«  ^^^^  ^^«  '-''^  ■  ^"■-'^ 

This  last  property  (Al-8)  makes  the  APW  (Al-l)  continuous 
at  the  sphere  radius  R. 

The  contribution  to  the  overlap  integral  ^k/k'^from  the 
volume  outside  the  spheres  is  just 


<^lt> 


and  that  from  inside  the  spheres  is 


<Klt>,-„   =  22  ''-««^^^  «/™f^';/«;(^,s)'-Vr  (ALIO) 

using  (Al-6).  This  can  be  simplified  by  applying  (Al-7)  to  derive 


the  relation 


2  <^*"^  *-^«.^^''  =  ^'''f^^^v^^*^^:^^*''^^  (^.11) 

".^ ^  ■  J'j,Ui')/u}iit.,e). 

Similarly,  for  the  energy  matrix 


(Al-13) 


because  of  the  relation  (Al-3). 


In  addition,  there  is  a  sui-face  contribution  due  to  the 
discontinuity  of  the  derivative  of  the  APW    |  kj>    at  r  =  R  ^ 

<  -t  /  <w  /  a:':>„^  =  -J  4^  8 1 74^,] .  d  t       (^.n) 
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where  6[VpJ   is  the  discontinuity  in  the  gradient  on  the  sphere  S. 
This  reduces  to: 

<tl  ^1 1>         -     ^^ZZ  ^-L^'^  ^^J^'^  "-p  ]  (Al-15) 

In  summary,  the  requisite  matrix  elements  between  two  single  APW's  are 

<tity  =  <k:ioiiKy  ^   Z"^^' ^■^'^>^^  (^-^^^ 
<:Si^iC>=  tt<k:i<ik.y 

^         _  -*  ^  (Al-17) 

■v^ere  Ji-  o 

~lZ-t'\  '        ' 

\ilYj,ll'>  '   4-n{l£^-l)  ix'^k.K)  jj^Ck/R.)  J>^U-k) 

<A1-19) 


=      -    2.     Die)  (Al-20) 


j;uj=  o,(.).£y€^  =  e,c.;-f34^.(^;  ^^_^^^ 


^^(^^-    ^Yf.-^"^  "^^^'"^J 


^'"'^  (Al-22) 
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The  secular  equation  is  obtained  by  forming    < kjX-E/ k/^ 
for  all  k  and  k  having  the  same  reduced  vector  k^,  and  finding  the 
zeroes  of  the  resulting  determinant: 

1^   -    E(p(     -f.^^   D^(e)   jfj,|     =    0   ,  (Ai-23) 

where       <?/ja/^'>=      tZ'    <t\«lt>  , 

The  symmetrized  combinations  of  AFW  functions  corresponding 
to  the  p^l^  representation  of  ^^^J,  the  group  of  the  reduced  wave  vector 


k^,  are  defined  by: 


where  ')>t*/jj  is  the  matrix  element  of  the  group  element^  corres- 
ponding to  the  I^^  row  and  the  J^^  column. 

The  matrix  element  between  2  symmetrized  APW's  corresponding 
to  the  same  reduced  vector,  of  any  operator  0  idiich  commutes  with  all 
the  group  elements  can  be  shown,  setting  k  =  k^,  +  g,  k''=  t.^  +  "g^  to 
reduce  to 

Tip       ^^        "  (Al-25) 

where  G^  is  the  order  of  ^Ck,),  n^  is  the  dimension  of  the  representation. 

Thus  for  a  particular  value  of  P  and  I,  the  secular  equation 
will  have  elements  of  the  following  form: 

<t  p  IT  I  ^-£  I  S'  P  I  ry 

=  ^  ^  r;^^^.  [a-^t'-  £)  <^i^ i^^>  ^"""''^ 

JL-=o 
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We  note  that  this  expression  is  independent  of  the  row  subscript  I, 
which  implies  an  n  -  fold  degeneracy  for  each  energy  eigenvalues 
of  the  secular  equation: 


I  IB  -  E  A\  +  ^  D^(E)  $^  I  =  0  ,  (Al-27) 

Jt 

P 


where 


JL 


\r   =   ±^<  itly,|tPJJ^>    . 


•Kp 


APPENDIX  II 


RECIPE  FOR  THE  CONSTRUCTION  OF  A  LINEARLY  INDEPENDENT  BASIS  SET 

The  problem  is  to  choose  from  among  the  functions  ft  |k_  +  g^ 
a  linearly  independent  set  of  functions  for  a  particvilar  row  I  of  the 
representation  p  of  ^(£j ,  the  group  of  the  wave  vector  t^, 

IT*   "♦v  . 

I  ^  +  g/  IS  any  Bloch  function,  such  as  an  APW,  corresponding  to 

the  reciprocal  vector  Icq  +  g. 

The  projection  operators  0^.  as  defined  in  equation  (2,29) 
satisfy  the  orthogonality  relation 

^JL  ^m      ""       ^"  hc(   ^LM  ^fN   •  (A2-2) 

Consider  the  function. 

Since  g^  +*i  g  is  also  a  reciprocal  lattice  vector,  then 
O  *Jk+Q^is  a  member  of  the  basis  set. 

NOW,  fff,  ^,  I  ii.+  p  -  ^  P'^K\r  ^^  *^ ' VJ> 

J 


«■ 


2  r,(^.).,    C>l\^.^P    , 
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so  that  the  function  O'lj^jl'l/J)  is  linearly  dependent  on  the 
set  of  np  functions  (  Ol^^  l*^e^?^   ,  L  =  1,  2  ..  .Hp).  Therefore, 
the  first  ingredient  of  the  recipe  isj 

RULE  1:  If  the  function  |ko+  g^  is  chosen  to  be  in  the  basis  set,  then 
no  functions  of  the  form  ^.Ik^  +  "g^  are  to  be  chosen. 

This  takes  care  of  the  case  of  unequal  vectors.  Now  consider 
the  functions  vdth  fixed  g.  Let  tJt  =  {"^^  +-^y.   .  ^he  number  of 
linearly  independent  functions  O^jj  3^  ,  (j  =  1,  2  . .  .n  )  is  equal 
to  the  number  of  non-zero  eigenvalues  of  the  matrix 

"  (A2-4) 

=  _?  ^^  1  0^j,ljsy  (independent  of  I)  . 

lap 

Suppose  yr   is  invariant  under  a  subgroup  Of^     of  ^Ck,^) 

(  (n.^    is  the  stabilizing  group  of  '^   ),  i.e., 

If  d(.    E   ^^  ,  then  the  use  of  (A2-3)  gives 

SO  that  the  corresponding  matrix  satisfies  the  relation. 
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-^ 


IF       is  a  representa- 


where  ^  and  ^  /    are  elements  of  ^^ 

tion  for  ^(f^J,   and  therefore  is  also  a  representation  for  the  sub- 
group ^^  ,  vdiich  will  in  general  be  reducible.  Therefore,  there 
exists  a  transformation  S  which  will  reduce  Tp  to  block  form. 


$  Tp  S    =  A\(Rji)  = 


- 
(D 

< 

(D 

Let  A\  ~  5  A\S  be  the  transformed  matrix,  which  will  have  the 
same  eigenvalues  as  /\  ,  From  (A2-6),  it  follows  that: 


A  /X{.\)    -   A\C^y)AV 


(A2-7) 


r 

Consider  the  decomposition  of  A\  into  the  same  block  form  as  t^ 


A\'  = 


'^^^ 

A„ 

< 

A\„ 

A\^ 

< 

A,,,. 

Application  of  (A2-7)  gives  the  equations 

By   Schur's  lemma,   since  (A2-8)  holds  for  all  elements  ^ 

•Jf  m  '^  in'i 

and  ^.  belonging  to  Of^ ,   then,  either  A\^- =  (D  or   iTq.  -  "o.  -  "^ 
(the  totally  symmetric  representation  of  Wyr  ) . 
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Thus,  the  number  of  non-zero  eigenvalues  of  A\  is  equal  to 
the  frequency  of  ITj   in  (Pp  .  This  can  be  expressed  in  mathematical 
form  to  give 

RULE  2:  The  number  of  projections  of  the  function  1^=1"^^  +  "g^that 
should  be  included  in  the  basis  set  for  the  representation  Tp  is: 


where  ^  is  the  stabilizing  subgroup  for  ^  with  order  G  ^  ,  and 
Xp  (.^)  are  the  characters  of  the  p^^   representation  of  (3!(k  \ . 

At  first  glance,  this  theorem  would  seem  to  be  of  little 
use,  in  choosing  the  basis  set,  since  it  gives  on3y  the  number  of 
projections,  and  not  the  actual  columns  of  the  representation  matrix 
which  must  be  used. 

However,  due  to  the  inherent  symmetries  of  the  functions 
1^0  "^  ^y  >   SpC^jvery  often  vanishes. 

To  show  the  usefulness  of  (A2-9),  consider  the  following 
three  examples: 

(1)  '  p  =  "]^  the  fully  symmetric  representation  oi  Cn.Ck^},   then 
X^C^j)  =  1  '^^li   and  therefore  Si-^^h  1  for  all  the  functions,  "^    . 

From  which  it  follows  that,  every  function  |  k^  +  ~^^  must  be  chosen 
once  and  only  once  in  the  basis  set  for   f^  . 

(2)  Suppose  5^"  has  no  intrinsic  symmetry,  i.e.,^^^  consists  only  of 
the  identity  element  C^. ,  then: 
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Therefore,  a  non-synmetric  kQ+g  must  be  chosen  ru  times 

for    r'p   . 

(3)  Suppose  ir   has  the  full  symmetry  of  ^fC'*  i-e-j  ^i>  =  ^^^o) 

">-'       a^(4)-.  J-   ^x,*(^)   -    in 

Therefore,  a  fully  symmetric  kg+g  must  be  chosen  once  for 
the  representation  ' -j_  ,  but  not  for  any  other  representation.  This  is 
the  case  if  g  =  o,  and  Kq  is  a  point  inside  the  Brillouin  zone. 

In  svmimary,  we  have  the  following  recipe  for  constructing  the 
basis  set.  For  a  given  reduced  wave  vector  k  ,  this  is  to  be  repeated 
for  each  irreducible  representation,  "p  ,  of  ^  C  ^^  /  . 

1.  Form  a  list  of  vectors  k  +  g,,by  running  through  all  the 
values  of  the  reciprocal  lattice  vectors  g,  up  to  the  value  at  'v^ich 
the  series  is  to  be  truncated. 

2.  Starting  at  the  lowest  magnitude,  choose  a  vector  k  +g^, 

-*  -* 
cross  off  the  list  all  those  vectors  k+g  vrtiich  are  related  to  it  by 

a  relation  k^+g  =3^(lc+g^),  K,£^Ck,^)  axid   continue  in  this  fashion. 

3 .  Find  the  stabilizing  group  of  each  vector  w  left  on  the 
list,  and  from  the  character  table  of  ^(k.^) ,   find  the  frequency 
number     a^(if-)    ^      J-         ^     Xp^iSl)      ■ 

4.  If  a  ( )^  )  =  0  cross  ^  off  the  list.  If  a  (  ^  )  7^  0 
and  n^  >  1,  then  the  actual  projections  must  be  formed  to  determine 
the  linearly  independent  set. 

5.  As  a  final  check,  the  total  nvmiber  of  projections  of  all 
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the  representations  must  be  equal  to  the  original  number  of  vectors 
on  the  list. 

In  table  A2-1  is  found  a  list  of  the  ap  (  ^  )  for  the  points 
of  high  symmetry  in  the  Brillouin  zone  for  the  face  centered  cubic 

structure.  The  notation  used  for  the  representations  is  that  of 

39 
Bouckaert,  Smoluchowski  and  V/igner,  -^^   and  are  found  along  the  top 

of  each  table.  The  general  form  of  the  W  's  found  in  the  vector 
list  is  found  along  the  leftmost  colimin.  The  small  letters  along  the 
top  are  the  type  of  atomic  orbitals  compatible  with  each  representa- 
tion, e.g.,  in  the  first  table,  it  can  be  seen  that  a  d-orbital  is 
compatible  at  kg  =  (0,0,0)  with  the  two  representations  1^2.5"  ^^^  '\x 
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Table  A2-1.  Frequency  numbers  for  the  fee  structure 


Gajnma 

1 

2 

12 

25' 

15' 

2' 

1' 

12' 

15 

25 

(0,0,0) 

s 

d 

d 

f 

pf 

f 

0  0  0 

1  • 

0 

0 

0 

0 

0 

0 

0 

0 

0 

a  0  0 

1 

0 

1 

0 

0 

0 

0 

0 

1 

0 

a  a  0 

1 

0 

1 

1 

0 

0 

0 

0 

1 

1 

a  a  a 

1 

0 

0 

1 

0 

1 

0 

0 

1 

0 

a  b  0 

1 

1 

2 

1 

1 

0 

0 

0 

2 

2 

a  a  b 

1 

0 

1 

2 

1 

1 

0 

1 

2 

1 

a  b  c 

1 

1 

2 

3 

3 

1 

1 

2 

3 

3 

X 

1 

2 

3 

4 

1' 

2' 

3' 

4' 

5 

5' 

(0,2,0)*(ir/a) 

sd 

d 

d 

f 

f 

pf 

d 

Pf2 

0  a  0 

1 

0 

0 

0 

0 

0 

0 

1 

0 

0 

a  0  a 

1 

0 

1 

0 

0 

0 

0 

0 

0 

1 

a  b  0 

1 

1 

0 

0 

0 

0 

1 

1 

1 

1 

a  0  b 

1 

1 

1 

1 

0 

0 

0 

0 

0 

2 

aba 

1 

0 

1 

b 

0 

1 

0 

1 

1 

1 

a  b  c 

1 

1 

1 

1 

1 

1 

1 

1 

1 

2 
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Table  A2-1  (continued) 


L                       1 

2 

1' 

2' 

3 

3 

(l,l,l)*(n/a)     sd 

£ 

P^ 

d2 

pr 

a  a  a                 1 

0 

0 

1 

0 

0 

a  a  b                 1 

0 

0 

1 

1 

1 

a  b  c                  1 

1 

1 

1 

2 

2 

W  1  1'         2  2<         3 

(l,2,0)*(Tr/a)     sdf      d  f      pdf      pdf2 


abO  1  0  0  1  1 

aOb  1  0  0  1  1 

abc  1  1  1  1  2 


Lambda       12    3       Sigma        1    2    3    4 
(x,x,x)     spdf^    f  pd2f2      (3,^x^0)     gp^2^   ^^   p^^2  pdf2 


a  a  a 

1 

0 

0 

a  a  0 

1 

0 

0 

0 

aba 

1 

0 

1 

a  a  b 

1 

0 

1 

0 

abc 

1 

1 

2 

a  b  0 

1 

0 

0 

1 

abc 

1 

1 

1 

1 

B3 


Table  A2-1  (continued) 

Delta  1    2     2»    1'  5 

(0,x,0)  spdf    df    df  pdf2 

r                                                      OaO  10000 

aba  1          0           1          0  1 

Oab  1          1          0          0  1 

abO  1          1          0          0  1 

abc  1          1           1          1  2 


Z  12  3  4 

(x,2Tr/a,0)  spd^f^      df      pdf^    pdf^ 


a  b  0  10  10 

a  0  b  10  0  1 

abc  1111 


S  12  3  4 

(x,2Tr/a,x)  sd2f2     pdf    pS^fS      ^^ 


a  0  a  10  0  0 

aba  1010 

a  0  b  10  0  1 

abc  1111 
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APPENDIX  III 

APW  CHARGE  DENSITY 
The  electronic  wave  fxmction  corresponding  to  the  ith  row 
of  the  pth  irreducible  representation  of  the  group  of  the  reduced 
wave-vector  k  is 

vfliere  |k  PIJ^  are  symmetrized  APW's  as  defined  in  equation  (Al-24), 

pIT 
and  '^jT*-?  *^®  ^^®  elements  of  the  eigenvector  as  obtained  from  the 

secular  equation  (Al-26).  The  sum  is  carried  over  those  combinations 

of  lattice  vector  ^and  column  index  J  which  give  a  linearly  indepen- 

dent  set  of  symmetrized  ARV's,  as  described  in  Appendix  II,   N  r* 

is  the  normalization  integral,  which  will  be  derived  in  equation 

(A3-11). 

To  see  the  explicit  spatial  dependence  of  the  wave  function 
inside  the  spheres  we  can  express  (A3-1)  in  the  form 

The  coefficients  a*  are  defined  in  equation  (Al-2).  The  spherically- 


averaged  charge  density  inside  the  sphere  is 


=  i<^^,'>Z''^^j' 


(A3-4) 


ji  =  o  rvt--JL 
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where  we  have  used  the  orthogonality  of  the  spherical  harmonics  (Al-6), 

(A3-5) 
Using  the  definition  of  a^  in  equation  (Al-2)  and  the 
addition  theorem  (Al-7),  this  can  be  reduced  to 

^^t'  (A3-6) 

By  using  the  orthogonality  properties  of  the  matrix  elements,  the 
expression  inside  the  curly  brackets  becomes 

f  «, 

'    <K^p    Hj,    I  i^.M'  ^  ^^'>    '  (A3-7) 

where  the  ^^  matrix  is  exactly  the  same  quantity  that  was  calculated 
for  the  secular  equation  (Al-27)  . 

Therefore,  we  have  as  a  final  expression  for  the  inside 
charge 
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density  (spherically  averaged) 


J     V 


The  outside  wave  function  has  the  form 

^"^  ^t,  -  Z,  ^VJ  ^  :P  ^^^'^t  ^  .   (A3-9) 

The  average  charge  density  outside  the  sphere  involves  the  same 
calculation  as  in  the  derivation  of  (Al-9),  and  turns  out  to  be 

where  the  «<  matrix  was  also  calculated  for  the  secular  equation. 
The  nonnalization  integral  can  now  be  obtained  by  the  equation 

0  "^ 

0 

where    T.  (  E)  -  * 

0 

The  equations  (A3-8),  (A3-9)  and  (A3-11)  can  now  be  used  to 
find  the  "muffin-tin"  part  of  the  charge  density  corresponding  to  an 
APW  wave  function.  Once  the  secular  equation  (A1-.27)  has  been  diagonal- 
ized  to  find  the   ^  1^^^ J  »  we  have  only  to  sum  up  over  quantities 
that  are  already  available.  The  determination  of  these  quantities  for 


'^'"^^>.B)  J\nsE).^a, 
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a  suitable  number  of  occupied  levels  gives  a  total  charge  density 
from  vrfiich  a  new  one-electron  potentieil  can  be  determined. 

In  Section  3.2,  we  refer  to  the  amount  of  s-like,  p-like 
and  d-like  charge  inside  the  sphere.  These  terms  refer  to  the 
quantities, 

for  ^=  0,1  and  2  respectively.  The  total  charge  outside  the  spheres 
is  just 


APPENDIX  IV 

DERIVATION  OF  THE  COULOMB  POTENTIAL  FROM  A  GIVEN  CHARGE  DENSITY 

The  generation  of  an  averaged  potential  (muffin-tin  type) 

18 
from  an  averaged  charge  density  is  a  simple  electrostatics  problem 

The  equation  to  be  solved  is 

V\  =  -  8  irj>    ,  (^_^) 

which  is  Poisson's  equation  in  atomic  units,  iirfiere  f  =  fo  (constant) 
outside  the  spheres  and  p  =  p  (r)  (spherically  symmetric)  inside  each 
sphere,  subject  to  the  condition 

4ir  j   tr'j>(t)  cL^    ^  y>o-^c    "  ^  '       ,  (M-2) 

o 

R  is  the  sphere  radius  and  Jl    =  Jl  -    ^H      is  that  part  of  the  unit 

J 
cell  outside  the  sphere.  There  is  also  a  point  charge  -  Z  at  each 

lattice  site.  (An  electronic  chcirge  is  taken  as  positive).  Since 

(A4-1)  is  a  linear  equation,  we  can  express  the  solution  as  a  siun  of 

two  potentials 

V  -  V,  ^  v^ 

O-j^  consists  of  a  charge  -  Z+Q  at  each  lattice  site,  a  spherical  charge 
^(r)  - J>^    inside  each  sphere,  and  zero  charge  outside.  P2  consists  of 
a  charge  -  Q  at  each  lattice  site  and  a  constant  charge  Ji^  elsewhere 
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(cf.  figure  M-l).  Q  is  defined  so  that  the  total  charge  in  both 
problems  is  zero,  i.e., 


vdiere  C^  is  the  charge  outside  each  sphere. 


(A4-4) 


V 


Figure  A4-1  The  Evrald  problem 

We  now  consider  the  two  problems  separately: 

1)  the  spherically  symmetric  charge  j>-^  in  each  sphere  produces 
a  zero  field  outside  it,  since  the  total  charge  inside  is  zero.  There- 
fore V-j^  inside  each  sphere  has  contributions  only  from  the  charge 
density  inside  that  sphere.  Choosing  V^  =  0  outside  the  spheres,  we 
have  the  solution  for  V-i 


^1    -  ° 


■f  1 


t         '   (A4-5) 


where   (r(»-)=  4  TT  t-'"  j3  (  0  ,   0;(t-)  r  A-lTY'-f, 
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2)  the  charge  density  p  2  produces  a  non-spherically  symmetric 
potential,  vrfiose  derivation  is  less  straighf orward .  Formally  it  is 
given  by: 


V,(?)=  -  y     19.    +  covstant, 


-*  _>. 

■where  t  are  the  lattice  vectors.  However,  the  sum  over  t  diverges,  and 

the  constant  is  infinite.  To  avoid  these  difficulties,  we  can  use  a 
mathematical  trick  to  extract  a  convergent  series  from  the  two  infini- 
ties in  (M-6).  By  using  the  identity; 
)0> 


X  Tnr   J 


(A4-7) 
0 

we  can  express  the  sum  in  {AU-S)  as 

*  7/=i  J    <^«-     ^    ^  f-   constant.  (A4-8) 

The  lattice  sum  in  (A4-8)  is  periodic  in  r  and  can  therefore 
be  expressed  as  Fourier  series,  i.e.,  a  sum  of  reciprocal  lattice 
vectors,  g,  which  is  easily  shown  to  be: 


■^  ,-V--?/'^  ^Ay 


t 


,  _  e         ^     (A4-9) 
The  trick  now  is  to  notice  that  the  lattice  sum  on  the  left  hand  side 


of  (A4-9)  is  rapidly  convergent  for  large  u,  and  the  reciprocal  lattice 
s\am  on  the  right  hand  side  is  rapidly  convergent  for  small  u.  Therefore 
we  divide  the  integral  over  u  in  (A4-7)  into  two  parts,  introducing  an 


arbitrary  constant  £  : 


(A4-10) 
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The  infinity  in  this  expression  arises  from  only  one  term, 
i.e.,  that  corresponding  to  ;^  in  the  reciprocal  sum.  Subtracting 
this  term  off  leaves  a  finite  potential,  which  is  dependent  on  £  . 
However,  the  derivative  of  (A4-10)  with  respect  to  £  is  a  constant 
independent  of  r ;  ^  _.  t^  ,  i. 

^v^      .  fra   2'c-^'/''*''>"'"^  f  1^  ^  e-'  ''  ' 

''    "    -^^^/^^  ^  ?  (A4.11) 

where  we  have  used  the  relation  of  (A4-9).  Thus,  by  adding  the  constant 
2  TV  Q//i£*  to  V2,  we  will  have  a  final  potential  independent  of  £  j 

-^  j%     3^  f      ^'"^'      (M-12) 

where  erfc  is  the  complementary  error  function.  The  constant  £  is 
chosen  so  as  to  optimize  the  convergence  of  the  two  sums. 

Equation  (A4-12)'can  be  put  in  dimensionless  form,  by 
setting  g  =  ^K,  t  =  a?,  ?=  ax,  o<  =  £a,  where  a  is  the  lattice  para- 
meter  andi2=  a  3/4  for  the  fee  structure: 

^  iK-ll  '^^  (A4-13) 

t 

Bnpirically,  it  is  found  that  if  the  sums  in  (A4-12)  are 
truncated  at  the  6th  nearest  neighbors,  then  foro(=  3  the  error  in 
a  V2/2Q  is  less  than  10"^  for  all  values  of  r.  The  values  of  this 
function  for  the  fee  structure  are  tabulated  in  reference  18.  These 
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values  were  recalculated  during  the  course  of  this  work  and  were  found 
to  be  in  agreement  up  to  the  fifth  decimal  place.  V2  has  a  singularity 
at  the  origin  due  to  the  charge  -  Q  situated  there.  From  the  calcula- 
tion the  behavior  near  the  origin  vras  found  to  be, 

where  b  =  4.584862. 

The  spherical  average  of  the  potential  inside  a  sphere  (A4-13) 
will  be  just  the  expression  (A4-14)  plus  a  term  due  to  the  constant 
charge  density  p^ » 

^   '     ^  y         T        /-J  (A4-15) 

Since  (A4-12)  is  valid  as  an  expression  for  V,  for  all  values 

i'  ^ 

of  g  ,  it  holds  in  the  limit  as   £ — *"  ©* 

This  is  £01  exact  expression  which  does  not  diverge,  but  is  too  slowly 
convergent  for  actual  calculation.  However,  it  can  be  easily  seen  from 
this  form  that  the  average  value  of  V2  taken  over  a  unit  cell  is  zero. 
This  allows  us  to  calculate  the  average  of  Vg  in  the  region  outside  the 
spheres,  from  the  relation 

'0 


/: 


from  which  we  can  find 


KvJ"'- 


a  I  3a         z^*-*  J    a.  ^^  ^ 
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Substituting  into  (AA-17),  the  values  R  =  a/2V2  ,  Sl^   =  a^/4 
-  4w/3  and  the  value  of  b  from  (A4-14),  we  find  that  c  =  .8I6314. 

The  equations  (A4-15)  and  (A4-18)  give  us  an  averaged  expres- 
sion for  V2  of  the  desired  "muffin-tin"  form.  This  is  now  added  to 
(A4-5)  to  give  the  total  potential,  vrtiich  simplifies  to 


a 

irfiere  y-  ^ 

t-       ►"  J.      -t 


(A4-19) 
(kt 


0 


For  convenience,  we  define  as  our  final  potential,  putting  all  the 
constants  into  V°^^; 

where  c'=  ^ (b  -  c  -  TT)  =  2.41583,  and  Qq  is  the  charge  in  the  region 
outside  the  sphere . 


APPENDIX  V 

DETERMINATION  OF  THE  DENSITY  OF  STATES 
The  density  of  states  n(E)  corresponding  to  a  set  of  energy- 
bands  E(k)  is  approximately,  from  equation  (3.2) 


7i(E)    ^      —  y     ^('^^  (A5_i) 


^^  e1  £CJC)  ^  EfAE 

where  4  E  is  an  arbitrary  energy  interval.  Equation  (A5-1)  is  exact 
in  the  lijnit  as  ilE  ->  0.  W(k)  is  the  weight  assigned  to  a  point  if 
corresponding  to  the  number  of  points  equivalent  to  it  in  the 
Brillouin  zone. 

Since  E(k)  is  calculated  at  only  a  few  points  in  a  typical 
band  calculation,  some  type  of  interpolation  scheme  must  be  used  to 
find  E(k)  at  a  large  momber  of  points,  in  order  to  be  able  to  get  an 
accurate  expression  for  n(E)  in  equation  (A5-1), 

As  a  first  step  in  developing  an  interpolation  scheme,  we 
realize  that  E(k)  has  the  symmetry  properties  of  the  full  point  group. 
For  the  fee  structure,  this  is  the  full  cubic  group  0  .  This  means  that 
E(k)  can  be  expanded  in  terms  of  Kubic  harmonics  (table  A5-1) 

vhere  J-(k)   is  the  spherical  average  of  E(k),  and  a.j   are  constants  which 
can  be  used  to  fit  (A$-2)  to  the  calculated  points. 
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*■*• 
^ 


Table  A5-1,  A  list  of  Kubic  harmonics 


Definition;  ^  *^  ^  *"' 

=  ^ 

't 

^[(rf^ 

(if] 

7  =  2. 

Coefficients; 

Ao 

A2 

A3 

\ 

h 

h 

-£  =  4 

-1 
5 

1 

- 

- 

- 

- 

6 

22 

11 

n 

1 

- 

- 

- 

8 

39 

1^ 
13 

-28 
15 

1 

- 

- 

10 

2^2 
3553 

-210 
323 

6^0 
323 

19 

1 

- 

12 

96577 
2kl5 

2475 
IhliS 

-660 
437 

495 
161 

-66 
23 

1 

^  The  Kubic  harmonics,  as  defined  here,  are  unnonnalized,  but  are 
orthogonal  to  each  other  and  have  a  zero  spherical  average. 
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However,  in  a  case  vdiere  d-bands  are  involved,  there  is  a 
great  deal  of  crossing  of  the  bands  along  symmetry  directions,  making 
it  difficult  to  define  a  function  E(k)  in  a  single  band.  It  is  more 
accurate  to  treat  E(k)  as  a  multivalued  function  -whose  values  are 


the  roots  of  a  pol^omial; 


2 


^,  (  O  £  '  = 


0 


(A5-3) 


where  N  is  the  number  of  bands,  and  the  k^  will  be  functions  of  the 
k-vector.  By  fitting  these  coefficients  A.  to  expressions  of  the  type 
(A5-2),  and  then  finding  the  roots  of  the  equation  (A5-3)  at  a  large 
ntunber  of  points,  an  accurate  expression  for  (A5-1)  can  be  found. 
In  the  actual  programs,  the  calculated  points  were  fitted  by  a  least 
squares  method,  subject  to  the  constraints  that  V,  E(k)  =  0  at  the 
points  X  and  L,  and  at  the  point  W  if  there  is  no  degeneracy.  This 
property  is  automatically  satisfied  at  the  point  T  because  the  Xji^ 
are  even  functions  of  r. 

The  calculation  of  the  density  of  states  curves  was  tried 
for  three  different  sets  of  k-points: 

1.  grid  #  3,  spacing  =  T/4a,  2048  points. 

2.  grid  H   4,  spacing  =  TT/Sa,  16384  points. 

3.  grid  #  5,  spacing  =Tr/l6a,  131072  points. 

It  was  found  that  it  was  necessary  to  use  the  largest  mesh  in  order 
to  reproduce  the  proper  behavior  around  the  van  Hove  singularities,^ 


APPENDIX  VI 

NOTATION 
a  lattice  parameter 

f{  sphere  radius  =  a  /  2  -/J? 

p  charge  density  for  spin  s  in  atom 

S  charge  density  for  spin  s  in  solid 

<'5j>        -spherically-averaged  density  in  solid 
V^   Vf.   <C^^>    Coiaomb  potentials 
\s  \s  <%s>    exchange   " 
\.^  intersphere  potential 

i.  lattice  vectors 

c 

7t.  number  of  neighbors  at  distance  a. 

Q  charge  of  spin  s  outside  spheres 

J2  unit  cell  volmne 

J 
J2  voltmie  outside  spheres  =  i2  -  f-F/y 

3 
H  ■  one  electron  orbital 

p.  radial  part  of  orbital 

g,,  one  electron  energy 

yi  Jl  S  atomic  quantiom  numbers 

<^( ''-'  radial  charge  density  =   4Tr>-^p(»-) 

^>iJLs  degeneracy  of  atomic  orbital 

j-  fraction  of  electrons  with  spin  s 
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jjy.  one  electron  wave  function 


Ej^  one  electron  energy 

K,  reciprocal  vector 

3  reciprocal  lattice  vectors 

K,.  reduced  wave  vector 
o 

Ce  -¥  components  of  eigenvector 

^_.j/  energy  matrix  element  between  basis  fvmctions 

A  J  J/  overlap  matrix  element  between  basis  functions 

I  j^  y  basis  function  (or  APW)  corresponding  to  k 

a,  (.II)  coefficients  in  spherical  expansion 

UaC*"   ^)  radial  function 

y  spherical  harmonic 
jc  yi^ 

01  L^  )  group  of  reduced  vector 

^  group  element 

O.J  projection  operator  for  pth  representation 

r  (.^)  .  J.  matrix  element  - 

Ik    PI  jy  projected  basis  function  (AFW) 

Q.  order  of  01  (^k.^) 

Tip  dimension  of  P 

£Lx)  J=    /     X>  o 

T  *    0    X  $  o  step  function 

Pj^C^.B)  =rUj,(r,E) 

cJ^  (  ^y)  Legendre  pol^omial 

^'.    (x)  spherical  Bessel  function 

>^  sphere  radius 

<^k.l  ft)>    ±.      overlap  integral  outside  the  spheres 

Ol/l 


i. 
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<A./  't/>/„  overlap  integral  inside  the  spheres 

°^}J    ^}fji  symmetry  dependent  matrices  found  in  secular  equation. 

^JL  ^  ^)  logarithmic  derivative  of  radial  function 

Ij^C^)     C  integrals 

A     B    S ,  symmetry  dependent  matrices  found  in  projected  secular 

i      J  -*  equation 

'^/>  characters  of  pth  representation 

2^   _i  one  electron  wave  function  corresponding  to  pth  repre- 

%  sentation  of  the  group  of  k 
f>fj  ° 

^  h.  +t  component  of  eigenvector  corresponding  to  basis  function 

(^Jy^  coefficients  in  spherical  expansion  of  the  wave  function 


^!n  spherical  charge  density  inside  the  spheres 

4^^.  ^oot  average  charge  density  outside  the  spheres 

I  ti  "^i  linear  independent  set  of  basis  functions  I  k  +e  PIJ^* 

'''^^  "^  for  P  and  I  fixed  °     ''^ 

/  2-  two  potentials  in  Ewald  problem 

Q  charge  defined  in  Ewald  problem  -  f^Sl 

J^o  constant  charge  density  outside 

Jl  cell  volume  =  a^/A. 

a^  total  charge  density  outside  =  fo-^o 

f  Ewald  parameter 

L  0^  =  ag 

^  =  4.584862  ] 

C  =  .816314  y  constants  in 

C'  =  2.415831  =  -(2  [ly.c^-fjj  I      Ewald  problem 
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